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TECHNICAL MEMORANDUM X-53850 


OPTICAL PROBING OF SUPERSONIC AERODYNAMIC TURBULENCE 
WITH STATISTICAL CORRELATION 

PHASE I: FEASIBILITY 

SUMMARY 


The theory and qualitative experimental results assessing the feasi- 
bility of measuring the statistical properties of supersonic turbulence 
by statistically correlating signals retrieved remotely with a laser 
quasi-schlieren system are presented. Also, some data retrieved with a 
laser shadow-correlation system are shown. Cross-correlograms and auto- 
correlograms computed on-line show that both of these systems can be used 
to retrieve flow-related signals sufficient for computing accurate and 
reproducible "peaks 11 of correlation. Additive tracers were not intro- 
duced because the schlieren and shadow sensing modes were used. 

A statistical method for obtaining "one-shot" measurements of the 
decay history of turbulent structures in a stationary frame of refer- 
ence is introduced, and results of practical applications of two types 
of these techniques are shown. The "one-shot" methods represent the 
only means by which the turbulence decay history can be computed from 
the same statistical sample of data (i.e., the information can be com- 
puted from a single composite signal retrieved during one run of the 
facility). Theoretical analyses show that these "one-shot" techniques 
will also yield results in three-dimensional turbulent flow regimes 
with only minor modification of the beam arrangement (and in some cases 
with the insertion of time delays between signals) provided signal- to- 
noise ratios of the raw data time histories are not prohibitive. 

The application of parallel-beam geometry was used to increase the 
s ignal- to-noise ratio over that of crossed-beam geometry since only 
qualitative results were sought. Although the crossed-beam geometry 
should be used for quantitative measurements, the potential use of 
parallel beams for retrieving quantitative results from a restricted 
class of flows should be investigated. 

Electronically induced time delays are used as a means for (1) 
zoning the "one-shot" auto- and cross-correlograms, (2) avoiding "peak" 
overlapping, and (3) identifying the "peak." 



With respect to the objectives of this investigation, the quali- 
tative results presented in this report clearly indicate that optical 
remote probing of supersonic aerodynamic turbulent flows is feasible, 
without the use of tracers (using statistical correlation). Further, 
correlation "peaks 11 computed from signals retrieved with either parallel- 
or crossed -beams during this investigation were flow-related, very repro- 
ducible and readily identifiable. 

No discussion of spectra is presented and no attempt is made to 
experimentally establish what flow properties are measured. These should 
be the objects of a systematic test directed toward obtaining quantitative 
results . 
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I. INTRODUCTION 


The majority of fluid flows encountered by the aerodynamicis t are 
turbulent. Because of the complex nature of turbulence, the approach in 
the subsonic regime has been through the application of statistical 
methods to data obtained by solid probes inserted into the fluid. How- 
ever, in supersonic and hypersonic flows, these probes adversely influ- 
ence the structure of the turbulence. One way to circumvent this 
difficulty would be to develop a reliable remote sensing tool for meas- 
uring the statistical properties of turbulence which would neither affect 
the flow field nor be adversely affected by it. The laser quasi-schlieren 
and laser shadow- correlation systems may, perhaps, represent the first 
generation of such a tool. 

The specific purpose of this publication is two-fold: (1) to pre- 

sent the theory and qualitative experimental results verifying the feasi- 
bility of optical remote sensing in supersonic turbulent flows employing 
the statistical correlation technique, and (2) to introduce and document 
the initial concept, theory, and reduction to practice of a statistical 
method which permits a "one-shot" measurement of the decay history of 
turbulent structures from one composite" time history of flow information. 

The remote sensing tool used for this investigation employs two 
laser beams of light (A = 6328A) which are retrieved by photodetectors 
after being influenced by the turbulent field. The time histories of 
these electrical ac-coupled (fluctuating about a zero mean) signals are 
amplified and filtered, and statistical correlation methods are used to 
retrieve the desired flow information (e.g. , speed profiles, eddy life- 
times, turbulent length scales, and spectra). 

The statistical "one-shot" autocorrelation technique constitutes 
methods whereby multiple time histories of random data are combined to 
form a single composite time history and an autocorrelogram is computed 
from the composite in such a way that a maximum amount of statistical 
flow information is retrieved with a minimum amount of time, equipment, 
and cost. The two basic types of "one-shot" correlations are referred to 
herein as the "one-shot" autocorrelation and the "one-shot" cross-corre- 
lation. Theoretical discussions, experimental results, and practical 
ramifications of these concepts are delineated in the main body of the 
report. 


The word "composite" is used in this report to imply a time history of 
random data composed of the algebraic sum of two or more random time 
histories of data. 
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The experimental objective of this work was to obtain accurate, 
reproducible, and readily identifiable correlation peaks with signals 
remotely retrieved from a supersonic turbulent flow which could be 
related to the most probable transit time of the disturbances. Previous 
attempts to obtain correlations in the supersonic regime have been 
hampered by (1) unknown anomalies in the data acquisition system (see 
Appendix C) , and (2) the influence of facility-induced noise upon- the 
optical system. 

The influence of facility- induced noise was reduced by utilization 
of the 7-inch Bisonic Wind Tunnel at MSFC, which operates at a very low 
noise level. This facility is described in detail in Appendix A. 

In addition to reducing the facility- induced noise, the power signal 
to- noise ratio of the raw data was increased approximately an order of 
magnitude by placing the laser beams parallel in the "two-dimensional” 
turbulent boundary layer on a thin-plate model. The laser beams were 
parallel to one another in the horizontal plane and normal to the flow 
direction, thereby increasing the correlated signals between the beams. 
This correlation technique is analyzed theoretically in section 3.1. 


Strong evidence that local information can be successfully retrieved 
from supersonic turbulent flows by using crossed beams is described in 
section 3.4.4. Also, the "cross-beam” method was used to investigate the 
effects the window boundary layers and interaction zones had upon the 
correlations computed from signals retrieved with parallel beams. It 
was found that, for the particular model design and thus the flow field, 
the boundary layers on the test section windows had no significant 
influence upon the correlograms . However, the interaction of the window 
and model boundary layers apparently dominated the measurements when the 
schlieren sensing mode was used in combination with parallel-beam geom- 
etry. Nevertheless, this does not affect the conclusions of this investi 
gation because the interaction zones near the windows are supersonic and 
turbulent. Also, the use of the laser shadow-correlation system reduced 
the contribution by the interaction zones approximately an order of rnagni 
tude. By crossing the beams in the vertical plane, these contributions 
to the correlograms can be avoided, as shown in section 3.4.4 and 3.4.5. 

The correlograms computed during this feasibility investigation 
represent our first encouraging measurements made in a supersonic turbu- 
lent flow with or without tracers and conclude the first phase of the 
wind tunnel "cross -beam" program. This has been conducted as an MSFC 
in-house research program with existing equipment and support. 
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II. BACKGROUND 


In November 1968, a test was initiated to determine the feasibility 
of retrieving signals by optical remote sensing of supersonic turbulent 
flows for obtaining accurate and reproducible flow-related statistical 
correlations. The major problem in attempting to retrieve signals by 
this remote sensing technique has been the low signal- to-noise ratio. 

One of the major contributors to noise has been the facility-induced 
(mechanical and acoustical) excitation of the optical system and even 
of the flow field itself. Since the larger, more advanced type of wind 
tunnel and air jet facilities being used in previous tests produced high 
noise levels, identifiable flow-related correlations could not be obtained 
in the supersonic regime. Because of the relatively small amount of 
noise produced by the Bisonic Wind Tunnel (BWT) facility at Marshall, 
this facility was selected for this investigation to alleviate the noise 
problem. 

The feasibility test in the BWT was planned in two parts, the objec- 
tive of part one being to isolate the optical system from facil ity- induced 
noise. The possibility for successful isolation looked very promising 
during the initial investigations. It is difficult to accurately esti- 
mate the s ignal- to-noise ratios obtained in previous facilities since 
facility-induced noise levels were always larger than flow-related 
signals, and also because it could not be determined that the correla- 
tions of signals were flow-related. 

In December 1968, a ratio of flow signal to noise of approximately 
ten was achieved in the BWT. A single laser beam was passed through 
the region of interaction of the shock wave and the free shear layer at 
position 1 perpendicular to the flow (figure 2-1) . The resulting signal 
is shown in figure 2-2B. The laser beam was moved to position 2 in the 
re-circulatory region of the base. The signal at position 2 is shown in 
figure 2-2A. All settings on instrumentation were the same for both 
runs. The experiment was repreated several times with the same results, 
thus implying that the increase in "s ignal" was flow-related.* The test 
objective of part one was achieved. 

The test objective of part two of the feasibility test was to obtain 
reproducible and readily identifiable correlation peaks related to the 
propagation of disturbances in the n two-dimensional” supersonic turbulent 


Position 1 in figure 2-1 was suspected to be one of considerable 
activity, because shadowgraphs of such shock wave shear layer inter- 
action regions indicate this. 


5 




6 



A 



Time Delay (Milliseconds) 



Time Delay (Milliseconds) 


FIG. 2-2 ON-LINE SIGNAL FROM SINGLE 
LASER BEAM PASSED THROUGH THE FLOW 
SHOWN IN FIGURE 2-1 

(A) TOP: BEAM AT POSITION 2 

(B) BOTTOM: BEAM AT POSITION 1 
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boundary layer on a thin plate. An attempt to cross-correlate the 
signals from two parallel laser beams in the BWT was made during a 
demonstration on January 22, 1969. The laser beams were separated by 
4.5 inches and approximately 1/8 inch above the surface of the plate 
(see figure 2-3). The signals from the two beams were amplified and 
cross-correlated in an analog correlator. The resulting cross-correla- 
tion, as a function of time delay, was displayed on a scope (figure 2-4A) . 
The maximum correlation occurred at approximately 260 microseconds, which 
corresponds to a propagation speed of about 1440 fps .* The free- stream 
speed was approximately 1660 fps . Because of the nature of the run, a 
great deal of care in locating the beams in the flow was not taken. 
Therefore, the measurement revealed only that the peak occurred approxi- 
mately at the expected time delay. 

Later investigations clearly showed that this correlation, though 
relatively weak compared to those computed later, was related to the 
propagation of turbulence in the boundary layer on the plate and in the 
interaction zone of the plate and window boundary layers. Further, a 
considerable amount of water vapor was present which extinguished the 
laser beams by scattering (schlieren effects were also present). Never- 
theless, this very crude measurement represented the first encouraging 
result from our attempts to retrieve signals from supersonic turbulent 
flows . 

In November 1968 during part one of the feasibility test, an attempt 
had been made to cross-correlate the signals from one laser beam which 
had been split with a beam splitter into two beams of equal power (a 
second laser was not available at that time). The purpose of this run 
was to determine if there were correlated harmonic signals present in 
the flow like those obtained during the preceding test in the 14- inch 
Trisonic Wind Tunnel (see Appendix B) . The two beams were passes through 
the free shear layer of the wedge shown in figure 2-1 . The correlation 
at zero time lag (see figure 2-4B) was due partially to the correlated 
laser noise. The large peak at about 100 microseconds gives a speed of 
1080 fps for the beam separation of 1.3 inches. The available supporting 
equipment and the test objective of part one would not allow the pursuit 
of this interesting result. Also, the facility- induced noise had not 
been fully investigated. Thus, no significant level of confidence could 
be put in the probability that the correlation or any portion of it was 


The peak of the correlogram will be used to determine the speed of the 
disturbances. This is not necessarily correct because of the dependence 
of the correlation function on both space and time. In this report, 
qualitative results are sought, and therefore this simplification should 
be justified. 
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GENERATED BY THE THIN PLATE MODEL 
INSTALLED IN THE 7X7° BWT WITH 
THE MACH 2 NOZZLE 


Cross-Correlation Cross-Correlation 
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Time Delay (Milliseconds) 



Time Delay (Milliseconds) 


FIG. 2-4. A(TOP): CROSS-CORRELOGRAM FROM 

TURBULENT BOUNDARY LAYER 
(M = 2.0) 

B(BOTTOM): CROSS-CORRELOGRAM IN FREE 
SHEAR LAYER OF 2-D WEDGE 
(M = 2.0) 
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actually flow-related. In retrospect, it seems reasonable to suspect 
that there was a significant flow-related contribution to the correlogram. 
The indirect influence of this particular run upon the test results is, 
perhaps, of interest and will be discussed later. 

After the cross -correlation was obtained in phase two (see fig- 
ure 2-4A), it was necessary to determine the origin of the signals 
from which it was computed. It was evident from the shape of the sub- 
sequent cross -correlograms of the individual beams (figure 2-5) that 
there was little or no periodicity present in the signals. Furthermore, 
because of the near-zero value of these cross-correlations at zero time 
delay and the absence of correlated periodic noise, it could be deduced 
that the correlated facility- induced noise was very small. Finally, the 
shape of the autocorrelograms indicated that the correlated signals were 
of the wide band variety [4] . Two additional conditions had to be met 
in order to establish the fact that the correlations were truly flow- 
related: (1) The peak had to occur at a time delay which corresponded 

to the approximate expected speed of the disturbances, 


U = — 

% 


and (2) condition (1) had to hold as the beam separation (|) was varied. 


The first runs to check these necessary conditions were made with 
the upstream beam (beam 1) and the downstream beam (beam 2) located in 
the wake of the model at positions 3 and 4, respectively (figure 2-3). 

The resulting cross -correlogram of the signals is shown in figure 2-6A, 
The beam separation and the time delay corresponding to maximum posi- 
tive correlation gave a speed of 1415 fps , which was at least reason- 
ably close to what was to be expected (free-stream speed was approximately 
1660 fps). Then, beam 2 was moved upstream such that the beam separa- 
tion (g) was decreased by one inch. The cross -correlogram for this case 
is shown in figure 2-6B. The maximum correlation corresponded to a 
speed of 1409 fps. Also, as would be expected, the strength of the peak 
correlation increased. The only difference between these runs, other 
than beam separation, was that the amplitude scale (voltage) on the 
oscilloscope was changed for the second run to accommodate the increase 
in amplitude of the correlation peak. 

These results, verified by later experimental data discussed in 
the following sections, provide reasonable experimental evidence support- 
ing the acceptance of the flow-related nature of the correlation peaks. 
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FIG. 2-5. A (TOP) : AUTOCORRELATION OF SIGNAL ONE 

B (BOTTOM): AUTOCORRELATION OF 
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Several new experimental and theoretical techniques were employed 
during this investigation which are, perhaps, worthy of note: (1) 

application of parallel beams for remote probing of a ” two-dimensional" 
turbulent flow, (2) the laser quas i-schlieren and laser shadow- correla- 
tion modes of retrieving signals from a supersonic turbulent flow, (3) 
the "one-shot" auto- and cross-correlation methods for statistical 
analyses of random time histories, and (4) the theory of induced time 
delays for peak identification purposes. 

Perhaps the most significant product of this feasibility test should 
be attributed to the attempt to correlate the random fluctuations of two 
beams of light originating from a common source. In itself, the splitting 
of a laser beam is, of course, common practice, but not previously applied 
using the cross-beam method. Nevertheless, the use of parallel beams from 
a common source leads to the initial M one-shot” autocorrelation concept 
described in section 3,2.1, and then to the "one-shot” cross-correlation, 
peak identification by the method of electronically induced time delay 
and a visual on-line display of the "one-shot” auto- and cross-correlo- 
grams . These methods are discussed in later sections. 
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III. DISCUSSION 


3.1 Remote Sensing by Optical Correlation 
with Parallel Laser Beams 

A theory of remote sensing of clear air turbulence using parallel 
laser beams (or other parallel light sources) in the supersonic regime 
is discussed, and experimental results are presented. 


3.1.1 Theoretical Description of A Laser Quas i- Schl ieren System 

In figure 3-1 the beams from two lasers, and L 2 , are directed 
through the two-dimensional turbulent boundary layer shown in figure 2-3. 
The flow is supersonic (M = 2.0), and the boundary layer is fully devel- 
oped along the thin plate. The beams, which are parallel to one another 
and perpendicular to the flow direction, are separated in the flow direc- 
tion (along the x-axis) in such a manner that the eddies which intersect 
beam LjlDi at time, t, also intersect beam L 2 D 2 at a later time, t + t. 

Since the flow is approximately two-dimensional, the time-averaged 
flow properties should be approximately equal across the test section 
(i.e., as measured along either laser beam). This does not necessarily 
require that the ins tantaneous flow properties across the test section 
be equal . 

The turbulence passing through one of these beams exhibits random 
fluctuating localized gradients of flow properties, e.g., local density 
gradient. As the eddies pass through a particular location on the beam, 
the resulting fluctuations of flow properties produce fluctuations of 
the local Poynting vector," s, of the laser beam. 

If j3(y , t) represents the local instantaneous angular deviation, 
as measured from a horizontal reference line perpendicular to the flow 
direction, then 


s 



cos p(y,t) > 


( 1 ) 


The Poynting vector is a vector tangent to the path of the beam with 
magnitude equal to the beam power. 

NOTE: The derivation presented here for parallel beams is similar to 

one presented for cross-beams by Dr. L. Wilson, IITRI, in a 
memorandum to Dr. F. R. Krause, MSFC, December 26, 1968. 
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FIG. 3-1. PARALLEL BEAM ARRANGEMENT FOR 
REMOTE ACQUISITION OF 7 INCH B.W.T. DATA 
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where j is a unit vector along the y-axis . Rearranging equation (1) 
results in 


P(y,t) = cos -1 



( 2 ) 


This fluctuation, j3(y , t) , is caused by a fluctuation of the local index 
of refraction, n (x, y , z , t) , and is related to the local fluctuation in 
the density gradient normal to the Poynting vector, s , for a particular 
wavelength of light. The local index of refraction, n(x,y,z,t), along 
the beam is proportional to the local density, p(x,y, z, t) . Provided the 
index of refraction is close to unity, as it is for most gases, the 
empirical Glads tone-Dale equation [1] is applicable. This equation, an 
empirical relationship between the index of refraction and the density 
of gas , is of the form 



Pi P2 


constant = a 


(3) 


so that 


n 0 (y,t) = l + a p Q (y,t) 


(4) 


for any point along the beam. 

Let us assume that n and p are represented by 


n(x,y,z,t) = n(x,y,z) + n' (x,y,z,t) (5) 

and 

p(x,y,z,t) = p(x,y,z) + p' (x,y,z,t), (6) 


where n' and p' are the fluctuations about the time-averaged values of 
n and p. Substituting equations (5) and (6) into (4) gives 


n = n + n' =1 +ap +ap' 


and 


= ap' 


(7) 

( 8 ) 
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It Is desired to relate the beam deflection angle, J3(y , t) , to 
the component of the index of refraction gradient which is perpendicular 
to the Poynting vector (path of beam). Since the index of refraction 
for air is very close to unity and since p(y, t) will be small. 


dp(y,t) 


1_ 5n 

n Bx 


dy. 


( 9 ) 


Equation (9), the derivation of which can be found In most textbook pre- 
sentations on the schlieren method of flow visualization [2] , represents 
the fundamental relationship upon which the schlieren method is based. 
Substituting equation (7) into (9) yields 


dP(y,t) + ^}dy. (10) 


Since n - 1 « 1 (n for air is approximately 1.00027), equation (10) can 
be simplified one step further: 

d ^ (y,t ) ^ lx dy + dy ' (11) 


The deflection of the Poynting vector away from the y-axis at a particu 
lar location on the beam (tj) will be obtained by integrating along the 
beam from y = - Z to y = q: 


P(r), t) 




+ 



Bn' 

Bx 


dy. 


-SL 


( 12 ) 


The first integral on the right-hand side of equation (12) represents 
the temporal -average angular beam deflection at q, and the second 
integral represents the fluctuation of the beam deflection about the 
temporal average. Thus, 


p(q,t) = p(ij) + -p' (q, t) . 


(13) 
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Substituting equation (7) into (12) provides a relationship between (3 
and p: 


P(t) > t) = 0! J dy + d; J dy. (14) 

-& -a 

The beam deflections in these relationships are in the xy-plane only. 

It will be shown later why the deflections in the z-direction can be 
neglected. Also, another assumption will be introduced: the magnitude 

of the time-averaged beam deflection, (3(i0, *- s such that the knife-edge 
(which is used to produce the fluctuations in current monitored by the 
photodiode) allows one-half of the laser beam light to pass into the 
photodiode when the beam deflection equals . Thus, equation (14) 

reduces to 



(15) 


This assumption is convenient and can be applied experimentally. The 
fluctuations of can be related to the fluctuating output signal, i(t) , 
of a photodiode, as shown in figure 3-1. 

A view taken along the laser beam (figure 3-2) shows the knife- 
edge, the eye of the photodiode, and the laser beam cross section. The 
beam is shown in a deflected position and in its time-averaged position 
(centered). Let 

t = time- averaged current output of the detector 
power supply when knife-edge is removed. 

A = cross-sectional area of laser beam. 

D = diameter of laser beam. 

A = beam displacement perpendicular to and measured 
from the knife-edge. 

AA = the change in area of the beam cut by a knife- 

edge due to a deflection of the beam off center. 
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Then, figure 3-2 shows that 


AA = 



(2A/D) + p 


VD 2 - 4a 2 \ 


(16) 


For small A, 


AA = DA, 


(17) 


The photodetector output per unit area of the laser beam is 


_d_ 

A 


fd 

JtD 2 


(18) 


To obtain the fluctuation of the output voltage (or current) 
from the photodetector, (i) , we now multiply the output per unit area 
by the decrease (or increase) in beam area cut by the knife-edge due 
to a beam displacement of A: 


i(t) 


4i,A 

d 

It'D 


for (A « D) . 


(19) 


Further, A is related to p 1 (J, t): 


A — d . 6’ (£,t). 


( 20 ) 


Substituting equation (20) into (19) gives 


i(t) = 


i-d I 

c 

JtD 


p' a,t). 


( 21 ) 
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The group of terms in parentheses on the right-hand side of equation 
(21) is a measure of the sensitivity, s, of the system [3], 


s 


4di d 
itD * 


( 22 ) 


The sensitivity of the photodiode with respect to beam position 
was assumed constant, as was the intensity profile across the laser 
beam diameter. These assumptions, however, may be far from true and should 
be checked. By placing a lens between the knife-edge and the photodiode, 
the beam movement on the photodiode can be considerably reduced. 

Substituting equation (22) into (21) gives 


i(x,t) = s • S' (x,i,t). 


(23) 


3.1.2 Cross-Correlation of Signals 

Equation (23) holds for both beams, so that 


= s 1 p' 1 .(x,,jg,t) 

i 2 (x + |,t) = s 2 • Pg(x + S,i,t) 


(24) 

(25) 


Substituting equation (15) into equations (24) and (25) provides rela 
tionships between i and cSp' /c)x, 


p '(x.ya.t) 

ii(x,t) = S ± a J 1 ^ " dy x 

-£ 


i 2 (x + |,t) = s^x 


5p 2 ( x + i»y2» t ) 

5x 


dy 2 , 


which hold along a particular eddy "streamline." 


(26) 


(27) 
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The signal, i 2 , from the downstream beam, is delayed by t seconds 
and the time-average magnitude of the cross product, l x • T^, is calcu- 
lated. This represents the cross-covariance of the signals I x and 
I 2 [I = 1/2 Tj + i(t) ] , or the cross-correlation of the fluctuating com- 
ponents, ii and i 2 . The cross -correlation of i x and i 2 will be used in 
the following discussion. 

The cross-correlation, R.(§,x), is defined as [4] 

T 

R(|,x) * lim 4 / ii(x,y 1} t) . i 2 (x + l,y 2 , t + x) dt, (28) 

T^oo L J 


where T is the averaging time. It is assumed that T is large enough so 
that 


R(£,t) = ™ / ii(x,y !, t) • i 2 (x + |,y 2 , t + x) dt 


(29) 


is a good assumption. Also, stationarity is assumed (i. e . , time-averaged 
statistical values are independent of time). 

In equation (29), R(g,x) will be maximum when i 2 is delayed by 
the most probable transit time of the common disturbances (x = x m ) pro- 
vided the signal-to-noise ratios and averaging time, T, are compatible 
and will allow the desired peak to rise well out of the correlated noise 
floor. Otherwise, the correlated noise will influence the shape and 
maximum position (x m ) of the peak. Correlated noise can often hide the 
peak altogether. The latter is likely to occur when the peak height is 
approximately equal to, or less than, the level of the correlated noise 
floor. 


Substituting equations (26) and (27) into (29) gives 


R(£,t) 



Sp^Cx, y ls t) 
Sx 




Spg(x + i,ys,t + t) 
Sx 


-dy ; 


t. 


( 30 ) 
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Rearranging the order of integration, we obtain 


R(|,t) = SjSzti? 


ax 



ap^(x+i 5 y 2 ,t+T) 

eSx 


dyidy 2 , 


■Z -Z 


(31) 


where the bar indicates the time-average of the product inside the braces. 
Thus, equation (31) provides one interpretation of what the cross-correla- 
tion represents. 

The cross-correlation of the signals i x and i 2 from two parallel 
beams is equal to a constant (s x s 2 a?) multiplied by the double line inte- 
gral, one taken along each beam, of the two-point-product time-average 
value of the fluctuating density gradient component in the flow direction. 

Another, and, perhaps a better, description is the one which 
shows the sequence of physical relationships between the parameters 
being measured experimentally. 

Beginning with the on-line analog correlator and moving opposite 
to the flow of data, the following physical events are taking place (see 
Appendix A, figure A- 11): 

(1) The scope displays the cross-correlation func- 
tion versus time delay (t) which is received 
from the analog correlator. 

(2) The correlator is computing the cross -correla- 
tion from the two electrical signals, i x and i 2 , 
originating from the photodiodes. 

(3) The time history of the signals are directly 
proportional to the angular deflections of the 
laser beams , respectively: 

i(t) = sp f U,t). 

(4) The angular deflection (£, t) of each beam is pro- 
portional to the fluctuating component of the 
density gradient in the direction of flow, 

where the mean is an instantaneous spatial 
average taken along the laser beam, i . e . , 
from equation (15): 
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(15) 


p' u, t) 


- <y. 

-i 


Multiplying and dividing the right-hand side 
of equation (15) by the width of the test 
section (2$ = L) gives 

- a d y} • 

-i 

Thus, at any particular instant (t ), 


S' <4,t p ) = 


5p (x,y, t p ) 


avg. over L 


(23) 


(32) 


The substitution of t = tp is made only to 
emphasize that the average is taken along the 
beam as if the flow were frozen in time. This 
is not implying that the beam length average 
value of SpV/Bx does not vary with time; it 
does. This results from the speed of light 
being much greater than the flow speed. Thus, 
the effect is (or can be) of considerable 
advantage to two-dimensional aerodynamic 
remote-sensing with statistical correlation. 

The advantage is a decrease in the integration 
time required to retrieve correlated statis- 
tical information from data in the presence 
of noise because of the automatic optical 
integration performed by the laser beam. For 
the two-dimensional case, the signals i x (t) 
and i 2 (t) represent the fluctuating com- 
ponent of the density gradient, (Sp* /5x), 
spatially averaged along the respective beams. 
Therefore, the correlator is operating on data 
which are representative of "instantaneous” 
averages of random samples of data. The con- 
vergence, therefore, will be faster than for 
the three-dimensional case where the beams must 
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be crossed. Further, and much more impor- 
tant, the signal -to-noise ratio is at 
least an order of magnitude greater for 
two-dimensional flow fields similar in 
flow rate and geometrical size. This is 
not meant to imply that optical remote 
sensing of three-dimensional flows is not 
possible. 

From the experimental point of view, these comments can be 
summarized as follows; 

(a) The correlation is representative of the similarity 
in the two electronic signals, i-j. and i 2 , received 
by the correlator after they have been influenced 
by the filters, amplifiers, and all other physical 
connections between the photo power supply output 
and the correlator input (see Appendix A, figure 
A-ll). If the correlator is performing properly, 
this is simply a matter of experimental fact. 

(b) The signals are proportional to the angular deflec- 
tions of the respective beams about a time-average 
value . 

(c) The signal represents the fluctuating (beam-average) 
x-component of the density gradient. 


3.2 The "One-Shot" Autocorrelation 
3.2.1 The Concept 

In figure 3-3, a single laser beam is passed normal to the flow 
and along the particular station of interest. After passing through 
the test section, the beam is reflected downstream from mirror Ml to 
mirror M2. From mirror M2 it is reflected back across the flow at a 
predetermined location and parallel to the first pass through the flow. 

A photodetector, Dl, receives the beam after its two passes through the 
test section. The signal is amplified and the analog correlator com- 
putes the autocorrelation versus time delay (autocorrelogram) , as shown 
at the bottom of figure 3-3. 

The autocorrelogram has two predominant peaks, one at zero time 
delay, as expected, and the second at a time delay, T m , corresponding 
to the most probable transit time of disturbances between the two posi- 
tions where the beam passes through the test section. Although the beam 
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M2 Mirror 
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arrangement is similar to that of the parallel beam case, the desired 
flow information is contained in a single signal rather than in two 
separate ones. The second peak occurs at x m on the autocorrelogram 
because the signal received by the photodetector, Dl, contains the 
same (or similar) information twice and at different times. The most 
probable separation in time corresponds to x m . 

This statistical technique can be used in analyzing random data 
such as those produced in turbulent fluid flows. The technique, 
referred to in this report as the "one-shot 11 autocorrelation, is 
described as follows: 

• The autocorrelogram of a single composite time 
history composed of the sum of two or more 
statistically correlated random signals which 
sufficiently lag one another will exhibit a 
correlation "peak" for each possible pair of 
the signals. 


3.2.2 Multiple Signals 

The extension of the concept would be to reflect the beam through 
the flow three or, perhaps, four times, as shown in figure 3-4. Although 
this technique will work in practice, it is preferable to avoid such a 
sensitive optical arrangement. 

In figure 3-5, an equivalent arrangement is shown. One laser 
beam is split into four beams of equal intensity and directed through 
the flow. Each beam has a separate detector that receives the signal. 

The signals are added and the resulting signal is equivalent to the 
original single beam case. It is not necessary to use a single light 
source. Each detector could have an independent light source, or two 
or more detectors may share light from the same source. 

The "one-shot" autocorrelation technique can, in theory, be 
applied to three-dimensional or two-dimensional flows. The only change 
required is a rearrangement of the beams so that no two are parallel 
and so that the effective correlated volume which is common to the beams 
is small enough to be acceptable. However, a longer integration time 
will be required for three-dimensional flows, as compared to two-dimen- 
sional flows, because of the lower signal- to- noise ratio of the former. 
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FIG. 3-4. MULTIPLE PASS PARALLEL BEAM ARRANGEMENT FOR 
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3.2.3 Peak Identification 


Consider a single time history of broad-band random data, F(t), 
(i.e., distributed over a wide range of frequencies) which is composed 
of two random time histories f x (t) and f 2 (t) such that 

F(t) = f x (t) + f 2 (t) (33) 

and the mean values of the time histories are zero. The autocorrelation 
of F(t) is 


T 

A(t) = lim ^ / [fi(t) + f 2 (t)]{f 1 (t + t) + f 2 (t + t)] dt (34) 

l-> 2 T-> oo 


A(t) = fx(t) • fx(t + T) + f 1 (t) • f 2 (t + T) + f 2 (t) • fx(t + T) 
1 ^ 2 


+ f 2 (t) • f 2 (t + t) (35) 


A(t) = A x (t) + R x 2 (t) + R 2 x (t) + A 2 (t) . (36) 

l-> 2 


In equation (36) letters "A" and "R" represent autocorrelation and cross- 
correlation, respectively. The subscripts indicate the signals involved 
in the particular correlation, as well as the order of the cross-correla- 
tion operation with regard to the time histories. 

Since 


F 2 l( T ) - Ri 2 (-t) > (37) 

Substituting equation (37) into (36) gives 

A(t) = A x (t) + R 12 (t) + R 12 (-t) + A 2 (t) . (38) 

1 -* 2 
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If we assume, for the present, that f 1 (t) and f 2 (t) are statistically 
independent and thus not correlated, equation (38) reduces to 


A(t) * A 1 (t) + A 2 (t). (39) 

1 ^ 2 


For this case, the "one -shot” autocorrelogram will have only one large 
peak* at t = 0 similar to the autocorrelogram shown in figure 3-6 (top). 
The magnitude of the peak at t = 0 is 


A(O') - A 1 (0) + A 2 (0). (40) 

1 ^ 2 


Beginning again with equation (38), we assume that f x (t) an d 
f 2 ( t ) are identically equal for all t when f 2 .(t) is delayed by r m sec- 
onds and compared to f x (t), i.e., at 


T 


m 


(41) 


f 1 (t) = f 2 (t m ) a for all t 


(42) 


where T m is the time lag of maximum correlation. For this case, 


F(t) - f x (t) + f 2 (t), (43) 

and substituting equation (42) into equation (35) gives 

A(%) = Ai(T m ) + Rl2(T m ) + R 12 (" T m ) + A 2 (T m ). (44) 

1 ^ 2 


The reference to the word "peak" implies the maximum positive peak in 
a local region on the autocorrelogram. The natural shape of the .auto- 
and cross -cor re lograms , near the time lag, r m , generally exhibits 
three or more peaks depending upon the frequency range of the corre- 
lated portion of the signals. 
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If we assume that T m is positive and large enough to prevent an inter- 
ference between the peak at T m and the zero time lag peak on the auto- 
correlogram, we can deduce from equation (44) that 


A (T m ) = Ri 2 (T m ),, (45) 

l-» 2 

since 

■ 0 

A.2'(%) = 0 

Rl2(-%) = 0 


(45a) 
(45b) 
(45 c) 


Equation (45) states that, if f-^t) and f 2 (t) are statistically random 
time histories and have maximum correlation when f 2 (t) is delayed in 
time by x m seconds, then the cross-correlation of two signals is equal 
to the "one-shot" autocorrelation provided that x is sufficiently 
large. Therefore, there should be no problem in identifying the peak 
at T m because there should be only one (see bottom of figure 3-6) . 
There is, however, a need to look a little deeper into the composition 
of the "one-shot" autocorrelogram. 

Consider the case of the addition of three random signals where 


F (t) = f x (t) + f 2 (t) + f 3 (t). (46) 


The "one-shot" autocorrelation of F(t) is 


A(x) - A^(t) + A 2 (t) + A 3 (t) + R xs ( t) + R 3.3 ( t) + R23( t ) ■*" R2l( T ) 
1-* 3 


+ R32(t) + R3i(t). (47) 


Using the general form of equation (37), we obtain 


A(t) = A 1 (t) + A 2 (t) + A 3 (t) + R 12 (t) + R 13 (t) + R 2 3 (t) 
1^ 3 
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+ RisC-t) + R£ 3 (-t) + R 13 (-t). 


(48) 



If it is assumed that all peaks resulting from the cross -correlations 
in equation (48) lie in the positive time lag range only, the last 
three cross -correlations in equation (48) drop out. 


A(t) = A 1 (t) + A 2 (t) + A 3 (t) + R 12 (t) + R 13 (t) + R 23 (t). (49) 

l-» 3 


Again, let us assume that the time lags between regions of correlation 
on the "one-shot" autocorrelogram are large enough to prevent overlapping 
of peaks. Let r 12 be equal to the time lag of maximum correlation 
between f x (t) and f 2 (t) , t 13 the time lag of peak correlation between 
f 1 (t) and f 3 (t ) , and t 23 the time lag of maximum correlation between 
f 2 ( t ) and f 3 (t) . Then, the value of the "one-shot" autocorrelation at 
the above will be 


A(t 12 ) — ^- 12 (^ 12 ) + ^23( t 12) 
L-» 3 

(50) 

4(Ti 3 ) = Rl3 (Ti 3 ) + R 2 3 (t 13 ) 
1-4 3 

(51) 

A(t 2 3 ) = R23 ( ^23 ) + Rl2 (^23) + ^13(^23) • 
l-» 3 

(52) 


If we assume that no overlapping of regions of correlation occur, then 


(50), 

(51), and (52) reduce, to 

A(t 12 ) 

= Rl2( T 12) 

l-> 3 


A(ti 3 ) 

= Rl3 (^ 13 ) 

l-s. 3 


A(T2 3 ) 

= R 23 (t 23 ). 

1-4 3 



(53) 


(54) 


(55) 
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Figure 3-7 shows a three-beam arrangement, where the parallel beams pass 
over the plate and through the turbulent boundary layer . The free- stream 
flow speed was approximately 1660 fps. The signals from the two photo- 
detectors A and B were added, amplified, and filtered. The resulting 
signal was autocorrelated on-line . The autocorrelogram is shown in 
figure 3-8 (top). The peak at the zero time lag is characteristic of 
all autocorrelograms of a wide-band time history. The peak at 
t = 0.128 ms is apparently a result of the correlation between beams 1 
and 2. It will be shown later that this peak is actually a result of 
the correlation between beams 1 and 2 summed with that between beams 2 
and 3. The peak at t = 0.254 ms represents the correlation between 
beams 1 and 3. The beam separations, (i.e . , | 12 is the distance 

between beams (1) and (2), etc.) were 

g 12 = 2.63" and g 13 = 4.95" 


which give average speeds of 


and 


5 12 


£ 

12t X 2 


1710 fps 


Ul3 


ii3 

12t 13 


1525 fps. 


The beams were approximately 0.2 inch above the plate, near the outer 
region of the boundary layer. The speed U 12 hardly seems reasonable 
because the free -stream speed was 1660 fps; the speed U 13 is closer to 
the expected value. However, there should have been a third peak repre- 
senting the correlation between the second and third beams, R. 23 ( T 23 ) • 
Dividing the beam separation, g 23 , by the speed, U 13 , should represent 
a good approximation to the location of the peak (t 23 ) (| 23 = i 13 - 
* 2.32"). 


g 23 

T 23 = = 0.127 ms. 

12013 

It appears that the "peak" at 0.128 ms is actually a result of the over- 
lapping of two peaks, and thus does not represent the independent loca- 
tion of either. To check the above logic, the first beam was moved 
1.00 inch downstream. This reduced the beam separation distances, 
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FIG. 3-7. PARALLEL BEAM ARRANGEMENT WITH 
THREE LASER BEAMS. BOTH, OPTICAL AND 
ELECTRICAL ADDITION OF SIGNALS ARE USED 
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FIG. 3-8.. "ONE-SHOT" AUTOCORRELATIONS AS 
OBTAINED WITH A THREE PARALLEL BEAM 
ARRANGEMENT AS SHOWN IN FIGURE 3-7 
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and !i 35 to 1.63 and 3.95 inches, respectively. However, the separation 
between the second and third beams, was not changed. The "one-shot" 
autocorrelogram for this beam geometry is shown in figure 3-8 (bottom). 

From figure 3-6B, we see 


t 12 = 0.082 ms 
t i3 = 0.203 ms 
t 23 = 0 . 1 23 ms « 


The corresponding speeds are 

U 12 = 1660 fps 
U 13 “ 1521 fps 
U 23 = 1570 fps. 


The purpose of this run was to show that the peak at r 23 did 
exist. The value of U 12 should not be considered accurate because the 
peak at t 12 = 0.123 ms is too close to r 12 to be considered accurate. 
There are, however, at least two ways by which peak overlapping can be 
avoided. These will be discussed in section 3.2.5. It is important 
that the number of peaks expected to appear for a particular number of 
beams is known in advance, as well as their expected location on the 
"one-shot" autocorrelogram. 

In general, the number of peaks on the "one-shot" autocorrelogram 

will be 


m-1 


No. of peaks = 1 + ) (m - k) 
k=l 


(56) 


provided that none are destroyed by overlapping. The peak at zero time 
lag has been included in equation (56). 

The "one-shot" autocorrelation could be used to study the character- 
istics of the decay of turbulent structures from data obtained from a 
single run of a wind tunnel or jet facility. The peaks resulting from 
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correlations between the first beam and each of the downstream beans 
give a picture of the "lifetime" of the average eddy in a stationary 
frame of reference. In addition, the peaks resulting from correlations 
between the second beam and the downstream beams represent the eddy 
lifetime a second, time, but beginning at a different location (i.e., 

| 2 further downstream). Likewise, the series of peak associated with 
the third beam and each beam downstream from it provides still another 
history of the decay* Therefore, the "one-shot" autocorrelogram is 
composed of at least one decay history (two beams) and, where more than 
two beams are used, there will be m - 1 of such histories, each history 
beginning a distance downstream from the last equal to the beam separa- 
tion (see figure 3-9). 

If the assumption is made that the decay history is independent 
of position over the distance between the first and last beams, we see 
that all of these time histories fall on the curve of the first time 
history, which begins with the first beam (see figure 3-10). It 
follows, for this case, that a peak resulting from correlation between 
two downstream beams provides the same information as the correlation 
between the first beam and an additional beam having the same beam 
separation as the two downstream beams. 

This characteristic of the "one-shot" autocorrelation is important 
with respect to the number of beams required to determine the decay 
history. For example, three beams provide four points on the decay 
curve, and four beams provide seven. 


3.2.4 Interpretation of the Zero Time Lag Peak 

The zero time lag value of the decay curve is not represented by 
the zero time lag value of the "one-shot" autocorrelation. However, 
this value can be calculated to a good approximation in many cases. 

Consider the general case of m number of time histories where 
f i (t) represents the ith time history and 

m 

F(t) =^T qct)- 

i=l 


The general form of the "one-shot” autocorrelation is 


T m 

A(t) = lim 4 / / f-(t) 

l-» m T-» oo 1 ^ £j 

O 1=1 


m 

^ f j (t + t) dt. 

3=1 


(58) 
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One-Shot" Autocorrelation 
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DECAY ENVELOPE USING THREE DETECTORS 




Equation (58) can be represented more conveniently in terms of auto' 
and cross-correlations, as was done previously. 


act) 

l-> m 



m 


I V T) - 

j=i 


( 59 ) 


For the applications considered here. 


R xy (t) = Ry X (-t) - 0, (60) 

for X > Y, and 

R xx (t) = A x (t), (61) 


where X and Y are particular values of i and j . Therefore, 


m 

A(t) = ) A. (t) .+ 

l-» m hi 
1=1 


m m 



i=l j=2 




j > t. 


(62) 


Equation (62) is the general form of the "one-shot" autocorrelation for 
this particular type of application. For the present, assume that the 
beams have equal sensitivities, s, and that the zero time lag value of 
the individual autocorrelations of the f^(t)'s are equal. Then, 


m 

A(0) = Y A i (0) = mA x (0) (63) 

i=l 

and 

A (0) = i [A (0) ] . (64) 

X m l-> m 
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If it is further assumed that the signal-to-noise ratio approaches 
infinity as the beam separation approaches zero between the first beam 
and any one of the downstream beams (i.e., no noise when = 0), then 
the value of the "one-shot" autocorrelation is equal to the cross- 
correlation at t = 0, 


A X (0) = R XY (T = °> ^XY = 0) * (65) 


For this particular case, the curve representing the trace of peaks as 
a function of t on the "one-shot" autocorrelogram (i.e., the curve repre- 
senting the rate of turbulence decay) will pass through the point 

[T = 0, i A(0) ] . 

1-4 m 


This value, at t = 0, is an upper limit and must decrease as the signal- 
to-noise ratio decreases, because, at t = 0 on any autocorrelogram, the 
noise is correlated as strong as the signal. 


3.2.5 Time Delays and Zoning of the n Qne-Shot fl Autocorrelogram 

Overlapping of regions of correlation on the "one-shot" auto- 
correlogram prohibits certain combinations of beam separations and 
limits the minimum separation between any two of the beams. There is 
a solution to the overlapping problem referred to herein as the method 
of induced time delay. Electronic time delay techniques are employed 
in many fields. In this report, there are three purposes of the induced 
time delay: (1) to avoid peak interference on the "one-shot" correlo- 

gram, (2) to identify a correlation peak partially lost in correlated 
noise, and (3) to zone the "one-shot" auto- or cross-correlogram, 
thereby providing means for identifying particular pairs of signals 
associated with a particular correlation peak on a "one-shot" auto- or 
cross-correlogram. 

With respect to (1) above, figure 3-11A shows a "one-shot" auto- 
correlogram where the measurement was made with two parallel laser beams 
separated by 1.90 inches in Mach 2.0 flow. The beams were located in the 
turbulent boundary layer. 

The effect of overlapping is clearly evident by comparing fig- 
ure 3-11A with 3-llB. In figure 3-llB, a cross-correlation was made 
between the two beams which eliminate the overlapping problem when only 
two beams were used. 
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There is a solution to the overlapping problem. If the second 
signal of figure 3-11A had been delayed by electronically introducing 
a known time delay, At, the peak at T m = 0.106 ms would have occurred 
at % = T m 4- At. If a time delay of 0.100 ms had been introduced into 
the signal from the downstream beam, the "one-shot" autocorrelogram 
would have looked like the one shown in figure 3-12, which is a much 
more accurate correlogram from the standpoint of studying the non- zero 
time delay correlation. 

In particular, induced time delays provide flexibility to the 
"one-shot” autocorrelation technique and extend its range of applica- 
tion, For example, the speed and direction of atmospheric ground winds 
can be theoretically determined by the "one-shot" autocorrelation method . 

Figure 3-13 is a schematic diagram of an atmospheric wind detec- 
tion system composed of six signals retrieved by optical remote sens ing 
techniques. Five of the six detectors are arranged in a fan for the 
purpose of determining which s ignal from the fan is correlated strongest 
with the signal from an "up- wind" location. 'It is also des ired to 
determine the time lag of maximum correlation. Ideally, this can be 
accomplished with the "one-shot" autocorrelation method using induced 
time delays of known value. 


As shown in figure 3-13, the signals from the fan B, C, D, E, and 
F are each delayed electronically and then added together. The result- 
ing signal is added to signal A. The "one-Shot" autocorrelogram of the 
summed signal should , theoretically, appear as shown at the bottom of 
figure 3-13. 

Induced time delays can be used to zone the "one-shot" auto- 
correlogram for the purpose of peak identification. If no time delays 
are introduced, the peaks would probably overlap. Even if no overlapping 
occurs, there would be doubt as to which signal from the fan produced a 
particular peak. 

The induced time delay, At, is large compared to the expected 
transit time , t m . Since signals B, C, D, E, and F are delayed by At, 

2 At, 3 At, 4At and 5 At seconds , respectively, a peak resulting from 
correlation between signals A and B would be expected to occur at 
T ab “ T m + At. Likewise, for signals A and C, t^c 58 % + ^At , etc. 

Also, the induced time delay method prov ides a graphical relationship 
between wind direction and the abscissa of the "one-shot" autocorrelo- 
gram. It follows that by comparing the "one-shot" autocorrelograms of 
the same data but with different induced time delays , peak identifica- 
tion should be enhanced. 

The induced time delay me thod of peak identification should apply 
just as well to cross -correlation, perhaps better . 
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3.2.6 Noise Analysis* 

The major disadvantage of the "one-shot" autocorrelation is that 
the noise-to-signal ratio, taken with respect to a particular time 
history contained in the composite, is greater than the noise-to-signal 
ratio of the particular time history taken separately. This is neces- 
sarily true, because, in a region on the autocorrelogram where two 
signals of a composite time history are correlated, all other signal 
and noise components represent noise in the computation. For applica- 
tions where the noise-to-signal ratios of the individual time histories 
are large, the "one-shot" integration time may become very long to 
achieve the same accuracy as that obtained from the individual cross- 
correlations. However, this has not been proven experimentally. 


In the following, let 

n a = noise-to-signal ratio of composite time history, 

n « noise-to-signal ratio of individual time history. 

m = number of time histories composing the composite 
time history. 

i = instantaneous value of the composite time history. 

i^ = instantaneous value of the kth individual time 
history (k * 1, 2, . . . , m) . 


i g = instantaneous signal of i^ 
k 

i n = instantaneous noise of i^. 


Assumptions ; 

(1) The noise-to-signal ratio of the individual time 
histories are equal. 

(2) The rms value of the individual time histories 
are equal. 

(3) There is statistical stationarity . 


The conventional signal -.to -noise ratio has been inverted here because 
the mathematics is simpler. 
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Derivation: 


From the definitions, it follows that 


m 


i 


c 


ii + i 2 + + . . . 



k=l 


( 66 ) 



= rms value of composite time history, 


[2 = 


JL 

T 


/ 


(ii + 


i ) ; 

m 


dt. 


( 67 ) 


i 2 

c 


i^ 4- i| 4* ... 4 2i^i 2 4 ... . 


( 68 ) 


The cross product terms approach zero as T a -> oo. Thus, 


m 

i 2 = if + if + . . . + ~\F = ) ip . 

c 12 m / , k 

k=l 


( 69 ) 


If the noise-to-signal ratio is defined as the ratio of the rms noise in 
the composite to the rms value of the signal, then, 


n 


a 



and 


( 70 ) 


n 




( 71 ) 
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If the individual noise- to-signal ratio is defined as the ratio of the 
rms noise to the rms signal, taken with respect to the individual time 
histories, then 


n 



and 



s 


Now substituting equation (69) into (71), we get 


(72) 


(73) 


n 

a 




By assumption 2, 



(74) 


(75) 

(76) 


and 





P + P. 

s n 


(77) 


By assumptions 1 and 2, 



(78) 



(79) 
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Substitution of equation (76) into (69) gives 



and from equation (77), we see that 


( 80 ) 



(81) 


Replacing ir in equation (80) by the expression in (81) yields 


i 2 = m(i 2 + i 2 ) . 
c s n ' 

Substitute equation (82) into (71). 



s 


Substitute equation (73) into (83). 

/■ -w 2 


n^ = -<jm(l + n 2 ) - 1 


(82) 


(83) 


(84) 


Equation (84) is the rms noise- to-s ignal ratio of a composite time 
history taken with respect to one of the individual signals and 
expressed as a function of the rms noise-to-s ignal ratio of the 
individual signal. 

The power noise/signal ratio of the composite time history 
is 


n = n 2 . 
pa a 


(84a) 
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The power noise/signal ratio of an individual time history is 


n = n 2 . 
P 


(84b) 


Then, from equation (84), we see that 


n pa = m(1 + n p^ " 1 * ( 84c ) 


In figure 3-14, n & is plotted versus m for several values of n. 
For large values of n, 


n = 'Jm' n. 
a 


(85) 


For the range of values plotted in figure 3-14, equation (85) is not a 
bad approximation even for n = 4, m = 2, For n = 0, which corresponds 
to the case where there is no noise in the individual time histories, 


n = 1 . (86) 

a 

It is seen from equation (86) that the minimum noise- to-signal ratio of 
the composite time history is unity. 

The practical limitation is not set by the increase in noise- to- 
signal ratio directly; rather, it is the increase in integration time 
required to obtain the same results that would be obtained by calculating 
the individual cross-correlations with equal confidence levels. The 
increase in integration time can be described as the ratio of T a to T. 
Here, T a is the integration time required to compute the "one-shot" 
autocorrelation and T is the integration time required to compute a 
cross-correlation between two of the individual time histories. The 
relationship is* 


h = 2 m 2 - 1 

T m “ (n 2 + l) 2 + 1 * 


(87) 


Vc 

This relationship is obtained from reference 6. The n used in the 
reference is the square of the n used here. 
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rms values And Noise - to- Signal Ratios of The 
Individual Time Histories Are Assumed Equal. 


Ratio of Noise- to - Signal of 
a Composite Time History, n a 


Number of Time Histories in Composite, m 


FIG. 3-14. NOISE -TO- SIGNAL RATIO OF A COMPOSITE 

VERSUS 

NUMBER OF TIME HISTORIES ADDED 
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In figure 3-15, the functional relationship between T a /T, m, and n 
is shown. From this figure the following observations can be made: 

(1) The minimum value of T a is 

T a = (2.5)1, 

which corresponds to the case of adding two identical 
signals to form the composite, i.e. , n = 0, m = 2. 
However, for this case, T is also a minimum. Therefore, 
adding two, three, or four such signals should be con- 
sidered practical. Verification of this has already 
been presented for m = 2 and m = 3. 

(2) The percentage of increase in integration time due to 
adding time histories with noise- to-signal ratios, n, 
equal to or greater than three (n § 3) is essentially 
the same for fixed m. That is, to a good approxima- 
tion, 

T 

-j = m 2 , for n s 3. 

(3) The limits T a /T are 


m 2 + 1 
2 


a 

T 


S m 2 . 


The increase in the required integration time, then, represents a very 
important limitation upon the practical application of the "one-shot" 
autocorrelation method. There is, however, another limitation that is 
not so critical. This is the increase in computation time. 

The ratio of computation times is given by 


7 


a (number of peaks ) • T 


( 88 ) 
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As explained in section 3.2.3, the number of peaks on the "one-shot” 
autocorrelogram is 


number of peaks 


m 


1 + 


(m - k). 


k=l 


( 56 ) 


The ratio of computation times is 


Thus , 


y* 




k=l 


(m - k) j 

-.1 


7 = 

■a 


m(m - 1) + 2 


T 

a 

T 


(89) 


(90) 


which has upper and lower limits: 


(m 2 + 1 ) ^ < 2m 2 

m(m - 1) +2 y a m(m - 1) +2 


(91) 


In figure 3-16, y a is plotted versus m. This figure represents the 
case where all peaks of the "one-shot" autocorrelogram are considered 
as useful data. 


3 . 3 The "One-Shot" Cross -Correlation 


3.3.1 Statement of the Concept 

The second kind of "one-shot" correlation technique is the "one 
shot" cross- correlation, which is described in the following statement: 
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n = co 



FIG. 3-16 RATIO OF COMPUTATION TIME ( / Q ) 

VERSUS 

NUMBER OF TIME HISTORIES IN COMPOSITE (m) 
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• The cross -correlation of a single time history 
with a composite time history, the single time 
history containing a random signal that is 
statistically correlated to two or more random 
signals contained in the composite, will result 
in a cross-correlogram exhibiting peaks of 
cross-correlation between the signal in the 
separate time history and each of the signals 
in the composite, the signals in the composite 
sufficiently lagging one another in time. 

The differences between the "one-shot" autocorrelation and the 
"one-shot" cross-correlation are as follows: (1) The "one-shot" cross- 

correlogram does not have the large peak at zero time delay representing 
the mean square value of the composite signal, and (2) the peaks result- 
ing from correlations between all possible pairs, other than those pairs 
formed between the signal in the separate time history and each of those 
in the composite, are not present on the "one-shot" cross-correlogram. 

These differences represent a loss of information. There are 
cases, however, where it is desirable to eliminate this additional com- 
putation. For example, when the noise- to -signal ratio of the individual 
time histories of a composite used in computing a "one-shot" autocorrela- 
tion are large, usually the error in evaluating the zero time delay point 
on the eddy decay curve is large. Also, the required integration 
time may be very long. In such a case, a considerable advantage can be 
gained by using the "one-shot" cross-correlation. The advantages are 
(1) lower noise -to --signal ratio and, therefore, a shorter required inte- 
gration time, and (2) the capability to compute the cross -correlation 
peak at zero time delay directly without an electrically induced time 
delay, which is necessary with the "one-shot" autocorrelation. 


3.3.2 Mathematical Development 

In the following, consider the cross-correlation of a single time 
history, f Q (t), with a composite time history, F(t) . F(t) is composed 
of m individual time histories, such that 


m 

F(t) =^T f k (t). 
k=l 


(92) 


Assume that each time history in F(t) and f Q (t) represents random time 
histories and that each can be expressed as the sum of signal plus noise 
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where a signal in F(t) represents that portion of a time history which 
is correlated to the signal of f 0 (t). Then, 


f D (t) - f So (t) + f„ o (t) 


( 93 ) 


m 

F(t) =y [f Si (t) + f n (t)i. 

Z_J k k 

k=l 


( 94 ) 


Furthermore, assume for each f^(t) that 


f k (t) 


J. 

- lim J f k (t) dt = 

T- — ) oo 


( 95 ) 


and 


T 


f o( fc > 


= lim — / f 

m T J o 

T— > <» 


(t) dt = 0. 


( 96 ) 


The "one-shot" cross-correlation of f Q (t) with F(t) is defined as 


R(t) «lim — / f (t) > F (t + t) dt 

k m T oo c ^ 
c o 


( 97 ) 


and 


m 


R(t) = lim f (t) 

l-> m T oo c ^ 
c o 


f k (t + t) dt - 


( 98 ) 


k=l 
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It is assumed that T c is sufficiently long, such that 


R(t) 
l-> m 



f k (t + T ) dt 


( 99 ) 


is a good approximation to equation (98) . 

Substitute equations (93) and (94) into (99). 

T 

c m 

r (t) = .£* / [f s (t) + f n (t)] . y [f (t + t) + f n (t + t)] dt. 

l->mc' JO o k 

o k=l 


( 100 ) 


Performing the indicated multiplication and taking the summation outside 
the integral, we obtain 


R(t) 
l-» m 



f s (t) • f (t + t) dt + / f s (t) 
o k d o 


f n, ( fc + T > dt 

k 


T 

c 


T 

c 


fn (t) 
o 


E a (t 
S k 


4- t ) d t + 


f n # f n & + T ) dt • 

ok 


( 101 ) 


The last three integrals on the right-hand side of equation (101) are 
approximately zero because the time-average value of the products of signals 
with noise, and noise with noise are zero; i. e . , the noise component of 
f 0 (t) is not correlated with the signal and noise components in the 
compos ite , F(t) . 

Equation (101) reduces to 


r(t) = y f f s ct> • f s (t + t) dt. 

l->m Cj i e « o k 

k=l o 


( 102 ) 
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The summation in equation (102) represents the sum of the cross-correla- 
tions of the single time history f Q (t) with each of the time histories 
in the composite evaluated at the same time delay, t. 


m 


R(T) 
l-» m 


= 1 E ok (T) - 
k=1 


(103) 


where the subscripts, ( ) Q | C > indicate the signals cross-correlated, as 
well as the order of the cross-correlation. 

If it is assumed that the signals in the composite are introduced 
sequentially such that each signal in time lags the previous signal suf- 
ficiently to avoid peak interference, the "one-shot" cross-correlation 
expressed in equation (103) becomes 


R(t .) = R ,(t ,), (104) 

o J oj oj N 7 

1— > m 


where t q j is the time delay corresponding to the peak of the cross -corre- 
lation Roj(T 0 j) between f 0 (t) and fj(t + t) in the composite. All other 
cross-correlations in equation (103) evaluated at T Q j are zero by the 
assumption that overlapping does not occur. 


I R ok<V ■ »• < 105 > 

k=l 


Therefore, the "one-shot 11 cross-correlogram will be composed of as many 
peaks as there are time histories in the composite, each of which con- 
tains a random signal statistically correlated to the random signal in 

number of peaks = m. (106) 

The time delay between the kth and (kth + 1) peaks is t q ~ t q 
which is sufficient to avoid overlapping. ? 
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3.3.3 Experimental Verification 

In figure 3-17, three laser beams, parallel to each other and 
normal to the mean flow direction, are directed through the turbulent 
boundary layer on a thin plate. The free- s tream Mach number is approxi- 
mately 2.0. The beams form a plane that differs only slightly from a 
plane parallel to the surface of the plate. Beam 1 from the upstream 
laser passes through the boundary layer, and a knife-edge blocks approxi- 
mately 50 percent of the laser light, preventing it from reaching the 
photodiode. The uninterrupted portion of the beam is monitored by the 
photodetector A. The fluctuations of the laser beam perpendicular to 
the knife-edge causes fluctuations in the amount of light monitored by 
the photodiode which converts these fluctuations into an ac-electr ical 
signal. 

The downstream laser beam is split into two beams of equal power. 

Beam 3 passes directly through the boundary layer. Beam 2 is reflected 
upstream and then through the boundary layer at a location which is 
2.32 inches upstream of beam 3 and 2.63 inches downstream from beam 1. 

After beam 2 passes through the boundary layer, it is reflected down- 
stream. Beams 2 and 3 are added optically by monitoring both beams with 
the same photodetector B. The fluctuations of beams 2 and 3 are retrieved by 
use of a knife-edge in the manner described from beam 1. 

The cross-correlation of the time history monitored by detector A, 
with the composite time history monitored by detector B, is shown in 
figure 3-7 (lower right). The two peaks predicted by the "one-shot" 
cross-correlation theory are clearly evident and provide strong experi- 
mental support to its validity. 

In order to compare the n one-shot n cross-correlogram obtained 
above with the "one-shot" autocorrelogram calculated from a composite 
of the same three time histories, the output of detector A was 
electronically added to the output of detector B, and the autocorrela- 
tion of the composite was computed as a function of time delay. This 
n one-shot ff autocorrelogram is shown in figure 3-17 (upper right). The 
structure of this "one-shot" autocorrelogram was previously discussed 
in section 3.2.3. 

In figure 3-18, a n one-shot ,f cross-correlation is shown with 
three beams in the composite time histories. During this particular 
run, the second beam in the composite was not positioned on the knife- 
edge. Instead, it was located on the edge of the photodiode opposite 
to the knife-edge (in the downstream direction). Thus, the edge of 
the photodiode performed the same function as the knife-edge, except 
the sign of the signal was reversed. This inverted the second peak on 
the "one-shot" cross-correlogram and demonstrates the schlieren principle 
very well. 
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MEASURED IN A SUPERSONIC TURULENT 
BOUNDARY LAYER 



0 ,1 .2 .3 .4 .5 


Time Delay (Milliseconds) 


FIG. 3-18. “ONE-SHOT" CROSS-CORRELOGR AM 
WITH THREE TIME HISTORIES IN THE 
COMPOSITE. SIGN OF SECOND SIGNAL WAS 

NEGATIVE 
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3.3.4 Noise Analysis 


For the same individual time histories the integration time 
required to compute the "one-shot" cross-correlation would be less 
than the integration time required to compute the "one-shot" auto- 
correlation, both having the same number of peaks on their correlograms 
(excluding the peak at zero time delay on the "one-shot" autocorrelo- 
gram) . 

In the following, we let 

= integration time for a "one-shot" cross-correlation. 

n = rms noise/ signal ratio of a signal in a composite 
time history. 

i (t) = instantaneous value of the independent time history. 

i (t') = instantaneous value of the composite time history. 

m » number of beams in the composite time history. 

1^.(1:) = instantaneous value of the kth individual time 
history in the composite. 

T = integration time required to compute the cross- 
correlation between the independent time history 
and one of the time histories in the composite. 

The rms noise/signal ratio of i Q is defined as 


n 

o 



(107) 


and from section 3.2, the noise- to -signal ratio for a signal in a com- 
posite is 


n c = [m(l + n 2 ) - l] 1 /- 2 . 


(108) 
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In equation (108), it has been assumed that the noise-to-signal ratio, 
n, and the mean squared values of the individual time histories in the 
composite are equal. Also, statistical stationarity is assumed. 

The ratio of integration times is 


T 

c 

T 


n 2 n 2 + n 2 + n 2 + 2 

c c 

n 4 + 2n 2 + 2 


(109) 


It has been assumed in equation (109) that the noise-to-signal ratio 
and mean squared value of the independent time history, i Q , are equal 
to those of an individual time history in the composite; i.e. , 


n = n 
o 


( 110 ) 


o 



( 111 ) 


Rearranging equation (109) yields 


T 

c 

T 


1 + 


(n 2 + 1) 

(n 2 + l) 2 + 1 


(n 2 - n 2 ) 


( 112 ) 


Substituting equation (108) into (112) and rearranging gives 


or 



1 + 


O 2 + 

-(n 2 +1) + 1 


(m - 1) 



■■ - 1? . 

(n 2 + l) 2 + 1 


(113) 


(114) 
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The upper limit of equation (114) is 


lim (T c /T) -> m (115) 

n-> 00 


and the lower limit is 


lira (T c /T) ^ (116) 

n-> 0 z 


or 


1 T 

m + 1 c 
2 - T 


§ m. 


(H7) 


Equation (114) shows that for n ^ 3 the ratio of integration 
times (T c /T) is, to a very good approximation, equal to the number of 
individual time histories in the composite; 


T 

"Y ^ m > for n ^ 3. (118) 

In figure 3-19, T c /T is plotted versus m for various values of n. 

From equation (56) in section 3.2.3, the number of peaks on a 
"one-shot 11 autocorrelogram are (excluding the peak at t = 0) 


m 

_ — j 


number of peaks = ) (m a - k) 

(119) 

L — ( 

k=l 


P a = \ (m a " m a>- 

(120) 


The number of peaks on the "one-shot" cross -correlogram are 



(121) 
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However, for comparison purposes, it should be required that 


P 

c 



( 122 ) 


Substituting equations (120) and (121) into (122) gives 


m 

c 



V- 


(123) 


Replacing m in equation (87) with m a , we obtain 


T 

a 

T 


m 2 - 
a 


m 2 - 1 
a 


(n 2 + l) 2 ■+ 1 


(124) 


which is the "one-shot" autocorrelation integration time ratio. Replac- 
ing m in equation (114) with m c , we obtain 


c 

T 


m 

c 


m - 1 

c 

(n 2 + 1) 2 ' + 1 


(125) 


Dividing equation (124) by (125) yields 


T m 2 (n 2 + l) 2 + 1 
a a 

T C m c (n 2 + l) 2 + 1 


(126) 


and substituting equation (123) into (126) yields 


T 2m 2 (n 2 + l) 2 + 2 

a a 

T c m a (m a - l)(n 2 + l) 2 + 2 * 


(127) 
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From equation (127), the following conclusions can be drawn- 


(1) The lower limit of the integration time ratio is 


lim (T a /T c ),-» 
n-> 0 


2m 2 + 2 

a 

m a (m a " 1 ) + 2 


(128) 


(2) The upper limit is 


2m 

lim (T a /T c ) -» — JY , (129) 

n-> oo a 


(3) The integration time required for the "one-shot" 

autocorrelation is always longer than that required 
for the ’’one-shot" cross-correlation. The minimum 
value occurs in the limit as m a -> oo and is equal 
to 2. 


T \ 
r). 

c-'mm 


lim (T a /T c ) -> 2. 

m — ^ oo 

a 


(4) The maximum value of T a /T corresponds to the limit 
as n for m a = 3 and is equal to 3. 


lim (T a /T c ) -> 3. 
n-> 0 

m a~ 3 



(5) For comparison purposes, the minimum value of m a is 
3 because the "one-shot" autocorrelogram for m a = 2 
has only one peak, other than that at t * 0, and the 
corresponding M one-shot" cross-correlation reduces 
to the standard cross-correlation. By equation (123), 
m c « 1 for m a = 2. These equalities imply that there 
is only one time history in the composite of the 
"one-shot" cross-correlation. ) 
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(6) The upper limit is approached rapidly as n increases. 
For n ^ 3, 


T 2m 

a a 


c a 

with less than one percent error. 
(7) The range of T a /T c is 

T 

2.50 g y- g 3.00 

C 


(130) 


for all combinations of n* s and m a f s between the 
limits defined below. 


0 g n ^ co 

3 ^ ^ oo 

cl 


Therefore, the required integration time for a "one-shot" autocorrelation 
is 250 to 300 percent longer than that required for a "one-shot 11 cross- 
correlation. This statement is, of course, restricted to the assumptions 
made above. 

In figure 3-20, T a /T c is plotted versus m a for several values of 
n. If the power noise- to-s ignal ratio had been used in the analysis in 
place of the rms noise- to- signal ratio, conclusion (3) would have been 
the only one affected (n i 3 would have been rip ^ 9 since n^ = n 2 ) . 

The ratio of the computation time required for the "one-shot" 
cross -correlation, T c , to that required for m number of separate cross - 
correlations (one for each peak on the "one-shot ,f correlogram) is 

7 = -£ = 1 - ' m . (131) 

C mT m[(n 2 + l) 2 + 1 ] 


The lower limit is 


lim (7 C ) 
n-^ 0 


m + 1 
2m * 


(132) 
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Number of Beams In Composite of 
"One-Shot" Autocorrelation, m a 


FIG. 3-20. "ONE-SHOT" AUTO- AND CROSS- 
CORRELATION INTEGRATION TIME COMPARISON 




and the upper limit is 


iim ( 7c ) -»1, (133) 

n— » oo 


so that 


na + 1 

~2S~ g 7 C 


s i. 


m ^ 2. 


(134) 


Thus, from equation (134), the computation time required to compute the 
11 one -shot" cross-correlation is always less than the time required to 
compute the m number of cross -correlations separately. However, as n 
increases, the ratio, y rapidly approaches its upper limit. For n ^ 3, 

7 c = 1.0, for n ^ 3 (135) 

with less than one percent error. Therefore, the "one-shot" computa- 
tion time is never a disadvantage. Actually, y Q does not account for 
computer printout time and other data handling procedures minimized by 
the "one-shot" method. 

Figure 3-21 shows y c plotted versus m for n = 0, 1 and oo. 

3 . 4 Cross Beam pi s cuss ion 


3.4.1 General Discussion 


The "cross-beam" theory [8] is essentially applicable to the 
measurements discussed below, which were made using crossed beams. 
However, as described in section 3.1.1, the signals were retrieved from 
the flow by employing the schlieren principle. The cross-beam runs dif- 
fer from the parallel beam runs in that the correlated signals from 
crossed beams are retrieved from a "localized" region defined by the 
segments of the two beams through which the same disturbances pass. The 
length of the beam segment is approximately the size of the average dis- 
turbance passing through and common to both beams. Because the beams 
are crossed, the flow disturbances which do not pass through both beams, 
and thus do not cause correlated fluctuations in the signals, represent 
noise. As a result, the cross-beam method should require more integra- 
tion time than the parallel beam method. Nevertheless, the crossed-beam 
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geometry offers many more advantages than disadvantages. The development 
of the cross-beam method is the ultimate objective* 

In the following discussion, a crossed-beam geometry was used 
to investigate the contribution to the correlations computed using 
parallel beams resulting from (1) the boundary layer on the test sec- 
tion windows, (2) the interaction of the wall boundary layer with the 
boundary layer on the plate, (3) the wake from the indentations in the 
plate running parallel to the windows (figure A-8) , and (4) a combina- 
tion of these. Also, a crossed-beam geometry was used to make measure- 
ments in the wake of the model. 


3.4.2 Effect of the Boundary Layer on the Test Section Windows 

In section 3 . 1.1, it was assumed that the statistical properties 
of the turbulence along each of the laser beams were constant. This 
assumption is not valid for positions near the test section windows 
because of the influence of the wall boundary layer. However, the 
magnitude of the influence was the important point to be determined. 

The necessity for this investigation results from the use of 
parallel beams. Because the average transit time of the disturbances 
between two parallel beams decreases from the window to the outer edge 
of the boundary layer, the effect should produce an unsymmetr ical con- 
tribution to the correlogram distributed from a time delay (t) slightly 
larger than the transit time corresponding to free-str_eam speed between 
the beams to larger values of t. Also, the total contribution from the 
wall boundary layer should be small compared with the correlation result- 
ing from the disturbances along the beams which are greater than the 
boundary layer thickness, away from the windows (see figure 3-22). 

Two runs were made using parallel beams separated by 1.87 inches. 
The first run was made with the beams passing through the boundary layer 
on the thin plate model (figure 2-3). The centerlines of the beams were 
0.05 inch above the surface of the plate. The cross-correlation of the 
two signals is shown in figure 3-23B. The time delay of maximum corre- 
lation, T m , is equal to 0.112 milliseconds (0,106 ms plus 0.006 ms due 
to parallax). The same result was obtained when the run was repeated. 

The beams were elevated to approximately 1.07 inches above the 
plate. One run was made and repeated twice. From the shadowgraph of 
the flow (figure 2-3), it can be seen that this is well above the 
major portion of the boundary layer (6 ^ 0.3 ,f ). 

All settings on the instrumentation were the same for all three 
runs made. The correlogram of one of the runs is shown in figure 3-23A. 
The correlograms of the other two runs were essentially the same as the 
firs t. 
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FIG. 3-22. SCHEMATIC OF FLOW IN TEST SECTION 

(PLAIN VIEW) 
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Time Delay (Milliseconds) 
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FI6. 3-23. A (TOP): CROSS - CORRELATION OF 

SEPARATED PARALLEL BEAMS, 
1.03 INCHES ABOVE PLATE. 
B(BOTTOM): CROSS - CORRELATION OF 

SEPARATED PARALLEL BEAMS, 
0.05 INCHES ABOVE PLATE. 
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A comparison of figures 3-23A and 3-23B would lead to the con- 
clusion that the boundary layer on the test section windows through 
which the parallel beams pass does not significantly affect the 
correlogram. The observation that the amplitude of the signals 
decreased by approximately a factor of five when the beams were in 
the elevated position well above the boundary layer supports this 
conclusion. 

Although the conclusion is reasonable, we cannot assume that the 
wall boundary layer 1.07 inches above the plate has the same effect on 
the correlogram as the flow resulting from the interaction of the wall 
and plate boundary layers near the surface of the plate. This region 
of interaction may contribute heavily to the correlation. Therefore, 
it was necessary to design a test which would separate the correlations 
along the beams, such that the contributions due to the interaction 
zone could be separated from those due to the boundary layer on the 
model . 


3.4.3 Investigation of the Boundary Layer Interaction Zone 

Figure 3-24 shows the plan view of the test section. The two 
laser beams were crossed such that the point of intersection was in 
the middle of the test section and the plane formed by the beams was 
approximately parallel to the surface of the plate. 

Two runs were made with this geometry and the cross-correlograms 
were computed on-line. Because of the geometry of the beams and the 
flow, the cross-correlogram is symmetrical about the origin (t = 0). 

The positive time lag range was computed during the first run (fig- 
ure 3-24A), and the negative range was computed during the second run 
(figure 3-24B). Cross-correlations between the beam segments AC and 
EC correspond to the positive range of the cross-correlogram, and cross- 
correlations between beam segments DC and EC correspond to the negative 
range. 


Two dominant peaks occurred on the cross-correlogram, one at 
+0.167 ms and the other at -0.167 ms (±0.158 ms ± 0.009 ms due to 
parallax). The common disturbances producing the correlations repre- 
sented by these large peaks can be estimated by multiplying the approxi- 
mate speed of transit, measured with parallel beams at the same distance 
above the plate (1510 fps) by the transit time (time delay) taken from 
the cross-correlogram. 


£ = U • T 

b approx. m 

= (1510) (12) (0.167) 10“ 3 
= 3.02 inches. 
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THE INTERACTION ZONE UPON PARALLEL BEAM CORRELATIONS 

A. (LEFT): POSITIVE TIME LAG RANGE 

B. (RIGHT): NEGATIVE TIME LAG RANGE 




Thus, the beam segments through which the common disturbances causing 
these large peaks pass are separated by approximately 3.0 inches. From 
the beam geometry, we see that the beam segments are very near the 
windows in the interaction zone of the wall and plate boundary layers. 

To check this result, the position of beam intersection was 
changed (see figure 3-25A). For this run, the beams were crossed in 
the wall boundary layer on the detector side of the tunnel. The plane 
formed by the beams was essentially the same as that of the previous 
case (0.20 inches above the surface of the plate). 

The cross-correlogram of beam 1 with beam 2 representing the 
cross-correlation between beam segments AU and CB is shown in fig- 
ure 3-25B. The time lag of maximum correlation is 0.263 ms 
(0.248 ms + 0.015 ms for parallax). Again using the mean speed of 
disturbances which was measured with parallel beams, the corresponding 
distance of separation can be estimated. 


S OQf . = U • T = (1510) (0.263) 1CT 3 (12) 
esc. m 

= 4.77 inches. 


Thus, the beam segments through which the common disturbances producing 
the correlation peak pass are separated by approximately 4.8 inches. 

From the geometry of figure 3-25A, we see that the beam segments 
separated by 4.8 inches, as measured in the direction of flow, are well 
inside the interaction zone near the window. 

The positive correlation at zero time lag in figure 3-25B was 
expected and is due to correlation of the disturbances passing through 
the intersection of the beams. 

Figure 3-26 shows the negative time lag range of the cross-correlo- 
gram representing cross-correlation of disturbances passing through beam 
segments DB and El. The zero time lag peak was expected; however, the 
large peak previously obtained in the negative time lag range with the 
beam intersection . in -the center of the test section was not computed. 

An acceptable explanation for this is not presently available. 

For the next two runs, the horizontal plane of the beams was 
elevated 1.55 inch above the surface of the plate. No other changes 
were made for these two runs, which were made to evaluate the con- 
tribution of the wall boundary layer to the correlogram. 
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FIG. 3-26. A ( TOP): SCHEMATIC OF BEAM GEOMETRY 
B (BOTTOM): NEGATIVE TIME LAG 
RANGE OF CROSS - CORRELOGRAM FOR BEAM 
GEOMETRY SHOWN IN SCHEMATIC 



Figure 3-27A shows the positive time lag range, and figure 3-27B 
shows the negative range. These cross-correlograms corroborate the 
results obtained using parallel beams, i . e . , that the wall boundary 
layer alone would have very little influence upon the correlograms com- 
puted from signals retrieved from the plate boundary layer if the two 
flows did not interact. 

Finally, it i can be concluded, but not without an element of 
uncertainty, that the interaction zone near the window and plate model 
is a region of intense turbulence and contributes significantly to the 
magnitude and shape of the correlograms computed from signals retrieved 
with parallel laser beams presented in this report. 

In general, these results do not invalidate measurements taken 
during this test series, since the test objectives were primarily to 
ascertain the feasibility of applying electro-optical remote sensing 
techniques to making measurements of supersonic turbulent flows, and to 
demonstrate the "one-shot" techniques without extensive alterations to 
the tunnel or expenditures of time and manpower. From this standpoint, 
the primary test objectives were achieved, since parallel beam correla- 
tion represents an optically integrated ensemble average correlation of 
"eddy” transit times across the entire test section. However, before 
quantitative measurements are made using parallel beams, a two-dimen- 
sional flow field with minimized edge effects should be sought, possibly 
through improvement in model design and sensing mode. 

Other attempts to minimize edge effects were also examined. One 
method investigated is a type of shadow-correlation technique, which is 
described in section 3.4.5. Crossing the beams in a vertical plane which 
is perpendicular to the flow also effectively eliminates edge effects, 
but the cross-beam method is rather cumbersome for measuring the velocity 
profile or other flow properties near flat plates because of geometrical 
limitations. Transparent models would, perhaps, allow cross-beam measure- 
ments, but passing a beam through a transparent model introduces a new 
problem of beam stability which should be investigated. Nevertheless, 
cross-beam measurements should be quite successful and convenient for 
making measurements in the wake. The next section discusses one such 
attempt to obtain these cross-beam measurements in the turbulent wake of 
the thin-plate model under investigation. 


3.4.4 Cross-Beam Measurement in Turbulent Wake 


This section discusses the feasibility of making cross-beam 
measurements in the turbulent wake of the thin-plate model. More space 
is available in the wake for geometric beam arrangements than in the 
boundary layer on the model. However, the beams could not be crossed 
perpendicular to one another [8] in a vertical plane perpendicular to 
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FIG. 3-27. A (TOP): POSITIVE TIME LAG RANGE OF 

CROSS-CORRELOGRAM FOR 

BEAM GEOMETRY SHOWN AT TOP OF 

FIG. 3-26. 

B(BOTTOM): NEGATIVE TIME LAG RANGE FOR 
SAME BEAM GEOMETRY AS IN (A). 
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the flow because the Bisonic Wind Tunnel does not have windows avail-* 
able for passing a vertical beam through the flow. Consequently , a 
modified cross-beam geometry was used (figure 3-28A). This beam 
arrangement was sufficient to avoid contributions to the cross-correlo- 
gram from the flow field near both windows. Only disturbances in the 
center of the wake and near the middle of the test section were common 
to both signals. The distance between the beams traversed by the common 
disturbances was 0.516 inch. 

The cross -cor re logram computed for this run is shown in fig- 
ure 3-28B . The maximum correlation corresponds to a time lag, T m , of 
0.036 ms (0.034 ms plus 0.002 ms compensation for parallax)* Thus, the 
most probable transit speed of a disturbance, averaged over the transit 
distance, is 


. TT . i 6 _ 0.043 „ 

< U> * rr TTv-g - 1194 fps. 

t 36 x 10 
m 


This speed is 71.9 percent of the free-s tream speed , U^, i.e. , 


< U> ^ 1194 

Uoo 1660 


0.719 


which is not unreasonable. 

The beams were not crossed in the vertical plane because it was 
desired that the distance between the beams vary across the test sec- 
tion. This provided means by which correlations near the windows, if 
they would have been present, could have easily been identified by 
their location on the cross -corre logram. 

This cross-beam measurement was apparently successful from a 
feasibility viewpoint. From these results, it appears that the appli- 
cation of the cross-beam method in measuring certain properties of 
supersonic turbulent flow is feasible. 


3.4.5 Experimental Results Using a Shadow-Correlation Method 


In the previous section, it was shown that the interaction zone 
near the test section windows had a significant influence upon the 
corre lograms computed from signals retrieved with parallel laser beams. 
The reasons for this are (1) the interaction zone is a zone of more 
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intense turbulence than is the turbulence in the plate boundary layer, 
and (2) the measurements were made with the laser quasi-schlieren system, 
which is sensitive to the first derivative of the density gradient com- 
ponent normal to the Poynting vector of the beams. 

There are at least four alternatives offering some relief to this 
problem: (a) to redesign the experiment such that the interaction zone 

is eliminated or the intensity of the turbulence in it is reduced to an 
acceptable level; (b) to employ a sensing mode which is less sensitive 
to the interaction zone; (c) to use crossed beams instead of parallel 
beams, or (d) to study the turbulence in the interaction zone itself, 
since it can already be measured very well. Because alternative (a) 
above is a matter of model or facility design, no discussion will be 
presented here. Alternative (c) has already been proven to be a 
successful solution (section 3.4.4). Alternative (d) provides interest- 
ing material for discussion but falls outside the scope of this report. 
Alternative (c) is of direct importance to the purpose here, but will be 
discussed briefly because of the limitations of the instrumentation 
available for investigation of the shadow-correlation remote sensing 
mode . 


In section 3.1.1, the laser quasi-schlieren principle was described. 
It was shown that the ac-coupled electrical time history of a beam was 
directly proportional to its deflection. This was accomplished by plac- 
ing a knife-edge in the beam (figure 3-1). As each disturbance inter- 
sects the beam, it causes the beam to be deflected in a particular 
direction, and the knife-edge is sensitive to the x-component of this 
deflection. Thus, it is evident that the detailed composition of a 
disturbance determines the magnitude and direction of a deflection. If 
this is actually the case, then it should be possible to sense the changes 
in composition by monitoring the intensity fluctuations inside the beam 
itself. To do this, the knife-edge in figure 3-1 was replaced with a 
thin plate having a small pinhole. The diameter of the pinhole was 
smaller than that of the laser beam. The plate is placed perpendicular 
to the beam such that the centerlines of the beam and pinhole are 
collinear (see figure 3-29) . Because the pinhole diameter is sufficiently 
small, the photodetector is never exposed to the circumference of the 
laser beam. Therefore, the fluctuations in total power monitored by the 
photodetector are a result of the fluctuations in intensity near the 
circumference of the pinhole produced by the disturbances inside the 
beam. 


Two runs were made using the shadow-correlation mode. The beam 
geometry for the first of these is shown in figure 3-30A. This run was 
made with the beams crossed in the hor izontal plane in the turbulent 
boundary layer on the plate model similar to those described in section 
3.4.3. The purpose of this run was to compute the contribution to the 
cross -cor re logram resulting from signals retrieved from the interaction 
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zones near the windows. The cross-eorrelogram is shown in figure 3-30B. 
At a time lag of 0.160 ms, the contribution of one interaction zone can 
be detected. However, by comparing this correlogram with that of fig- 
ure 3-24A, it can be seen that a significant reduction has been achieved. 

The purpose of the second run was to compute the correlogram from 
signals retrieved with parallel beams in the turbulent boundary layer 
for a qualitative comparison with figure 3-30B. The correlogram for this 
run is shown in figure 3-31B and the beam geometry is shown in figure 
3-31A. 


Comparing figures 3-30B and 3-31B, it can be concluded that the 
shadow-correlation mode offers a considerable advantage over the quasi- 
schlieren mode where the beams must pass through regions of turbulence 
near or on the test section windows. However, the magnitude of the 
correlation is reduced for the shadow-correlogram due primarily to a 
decrease in signal-to-noise ratio. 

These conclusions must be considered as preliminary due to the 
limitations of the instrumentation and hardware available. 
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IV. CONCLUSIONS 


Theoretical discussions and qualitative experimental results have 
been presented in an attempt to assess the feasibility of remotely 
retrieving statistically correlated signals from supersonic aerodynamic 
turbulence which could successfully be used to compute auto- and cross- 
correlograms exhibiting accurate, reproducible, and readily identifiable 
flow-related "peaks” of correlation. With respect to this objective, 
the qualitative results presented in this report indicate that optical 
remote probing of supersonic aerodynamic turbulence using statistical 
correlation is feasible, without the use of tracers, in the Bisonic Wind 
Tunnel of MSFC, and, perhaps, in other facilities which exhibit similar 
low facility-induced noise levels. Other conclusions of this investiga- 
tion are as follows: 

(1) Optical remote sensing with parallel or crossed laser beams 
can yield consistent and reproducible statistical results which are 
related to aerodynamic turbulence, the particular relationships requir- 
ing theoretical analyses and quantitative experimental verification. 

(2) The use of parallel laser beams for retrieving signals from 
two-dimensional turbulence* increases the power s ignal- to-noise ratio 
relative to that obtained with crossed beams. However, direct localized 
measurements require the crossed-beam geometry, 

(3) The laser quas i-s chlieren system presented herein requires 
isolation from the mechanical, and perhaps acoustical, facility- induced 
noise. For this system, the major portion of the facility-induced noise 
enters the system at the light sources. 

(4) With respect to the data presented in this report, which was 
obtained with a laser quas i-s chlieren system in combination with a paral- 
lel beam geometry, the boundary layer on the windows of the test section 
did not significantly contribute to the computed correlograms . How- 
ever, the interaction zone resulting from the interaction of the window 
and model boundary layers apparently dominated the measurements. This 
result does not invalidate measurements taken during this test, since 
the test objective was primarily to ascertain the feasibility of apply- 
ing electro-optical remote-sensing techniques to making qualitative 
measurements in a supersonic turbulent flow. From this standpoint, the 
primary test objectives were achieved. However, before quantitative 


Note: Two-dimensional time-averaged turbulent properties. 
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measurements are made using parallel beams, a two-dimensional flow field 
with minimized edge effects should be sought, possibly through improve- 
ment in model design and sensing mode, 

(5) The shadow-correlation method significantly reduces the con- 
tributions of the interaction zones , and the crossed-beam geometry 
successfully avoids them. 

(6) The "one-shot" auto- and cross -correlation methods provide a 
convenient means for measuring the decay history of turbulent structures 
in two- or three-dimensional fluid flows (for detailed comparison of 
these methods see section 3 , 3 . 4) . 

(7) The most important advantages of the "one-shot" autocorrelation 
method are as follows: 

(a) It represents the only method for measuring the decay 
history of turbulence from one sample of data. 

(b) Because it is necessary to monitor only a single 
composite signal, there is a reduction in instru- 
mentation compared to that required for a two- 
channel cross-correlation system. 

(c) Phase matching errors are eliminated. 

(d) Facility operation is reduced by 400 percent because 
only one run of the facility is required to measure 
the decay history. 

(e) The amount of data handling, storage, etc., is reduced. 

(f) The overall data-process ing time is reduced. 

(g) This method is the only means by which the decay history 
can be observed visually "on-line" by displaying the 
"one-shot" autocorrelograms on an oscilloscope as the 
data are collected and computed from the same set of 
turbulent data. 

(8) The major disadvantage of the "one-shot" autocorrelation is 
that the signal -to- noise ratio, taken with respect to a particular time 
history contained in the composite, is smaller than the s ignal-to-noise 
ratio of the particular time history taken separately. 
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(9) The "one-shot 11 cross-correlation method has the following 
major advantages over the "one-shot" autocorrelation: (a) There is a 

higher s ignal- to-noise ratio, and therefore the required integration 
time for a "one-shot" autocorrelation is 250 to 300 percent longer than 
that required for a "one-shot" cross-correlation. (This statement is, 
of course, restricted to the assumptions in section 3.3.4.) The "one- 
shot" method provides a means for computing the cross -correlation peak 
at zero time delay directly without the electrically induced time delay 
necessary with the "one-shot" autocorrelation. 

(10) The computation time required for the "one-shot" cross-correla- 
tion is equal to or less than that required to compute the equivalent 
number of individual cross-correlations. 

(11) The method of induced time delay extends the potential applica- 
tion of the "one-shot" correlation methods and may be used (a) to avoid 
peak interference on a "one-shot" correlogram, (b) to identify a correla- 
tion peak partially lost in correlated noise, and (c) to zone the "one- 
shot" auto- or cross-correlogram, thereby providing a means for identify- 
ing a particular pair of signals associated with a particular correlation 
peak on a "one-shot" auto-^ or cross-correlogram. 

(12) The laser quas i-s chi ieren system described in section 3.1.1 is 
the most reproducible, convenient, and flow-sensitive system presently 
available for studying remote sensing of aerodynamic turbulence in wind 
tunnels . 

(13) No experimental attempt was made to determine how the correla- 
tions computed from signals retrieved using the laser quas i-s chi ieren 
system were physically related to the flow disturbances causing the 
signals. It can be concluded, however, that the correlations presented 
in this report are related to the propagation of supersonic turbulence. 

(14) The shape of the correlogram in the region of a "peak" correla- 
tion computed from laser quasi-schlieren signals is directly related to 
the position of the knife-edges relative to one another and the flow 
itself. One explanation of this is presented in section 3.1.1 where it 
is theoretically concluded that the correlation between signals is pro- 
portional to the time-averaged two-point product of the fluctuating 
density gradient component in the direction of flow. For parallel beams, 
this two-point product is integrated along the beams from source to 
detector . 
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APPENDIX A 


DESCRIPTION OF THE BISONIC WIND TUNNEL 
FACILITY, MODELS, AND INSTRUMENTATION 


This appendix presents a general discussion of important design 
characteristics of MSFC 1 s 7-Inch Bisonic Wind Tunnel facility, the 
two-dimensional, boundary layer and wedge models, instrumentation for 
remote sensing supersonic flows, the general test procedure, and the 
analog data reduction equipment used for this test series. 

The salient features of this equipment and instrumentation and 
their respective roles in data interpretation are discussed below. 


A-l MSFC's 7 x 7 -Inch Bisonic Wind Tunnel (BWT) 

This test facility is a supersonic blowdown type of wind tunnel 
(Figure A-l). Dry air at atmospheric pressures (or atmospheric air at 
approximately 14.4 psia with 1.6 percent variation) is supplied to a 
7 ,r x 7 M test section and exhausted to essentially vacuum conditions 
(approximately 20 to 29 in. hg.). Figure A-2 shows a schematic repre- 
sentation of this facility. Major dimensions of the tunnel are given 
in Figure A- 3. 

Two air dryers are connected in tandem to charge the tanks (Fig- 
ure A-2). The first takes in air at 3,000 cu ft/min, while the second 
takes in this amount and recirculates an additional 12,000 cu ft/min. 

Two combined tanks equipped with a rubber diaphragm liner store up to 
60,000 cu ft; however, one tank is currently not in the circuit. 

Two Fuller duplex and two-stage rotary pumps provide an overall 
vacuum rate of 8920 cu ft/min. This results in a 40 -minute pump -down 
time from atmosphere (approx. 14.4 psi) to 0.15 psi for the six inter- 
connected cylindrical steel vacuum tanks. The combined capacity for 
this vacuum storage is 42,000 cu/ft. 

Mach number control, is provided by removable nozzles (Figure A-3) 
machined from solid brass to a tolerance of ±0.001 inches. A hydraulic- 
ally activated mechanism clamps and aligns the nozzle in position. The 
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test section formed by the nozzle is 7.029 inches wide and about 7 inches 
high (see Figure A-4) . Although several nozzles providing Mach numbers 
of 1.54, 1.99, 2.44, 3.00, 3.26, 4.00, and 5.00 are available, only a 
Mach number of 1.99 was used in these tests. The test rhombus, shock 
angle, mass flow, etc., for this nozzle-model combination are discussed 
in subsection A-2 of this appendix. 

The test section side walls which enclose the windows are hinged 
at the top and are activated by (see Figure A-l) hydraulic cylinders. 
These doors, one on either side, extend the full length of the test 
section. The 6- x 6 - inch glass windows, 1 and 1/16 inches thick, are 
mounted in the doors in such a manner that discontinuities between frame 
and glass surface inside the test section (Figure A-4) are avoided. 

These windows provide schlieren and shadowgraph visualization of the 
flow region, as well as to facilitate access for remote sensing by 
lasers, ultraviolet, or other light sources. 

Pneumatically inflated seals are used along the nozzle contour 
(Figure A-4) and at the rim of the doors. Tunnel operation is auto- 
matically blocked if all seals are not activated. 

Two primary valves (Figures A-2 and A-3) separate the test section 
from the air supply and vacuum tanks, respectively. Both valves, which 
are hydraulically activated, provide positive sealing when closed. When 
in the open position, a box-shaped extension on the gate is aligned with 
the tunnel duct in such a manner that a smooth continuity of aerodynamic 
surfaces is provided. 

Starting and stopping processes require about four seconds. The 
air flow is controlled by the upstream gate valve. This has two 
advantages: (1) starting loads on the model and test section are 

reduced to a minimum, and (2) the tunnel is evacuated before the flow 
is established, and therefore settling time is reduced for the flow 
since test section pressures are nearer vacuum than atmospheric. 

An adjustable outlet diffuser (Figure A-3) of high aerodynamic 
efficiency increases tunnel running time. The maximum and minimum 
positions of the diffuser are preset for each nozzle (corresponding to 
a given Mach number). At the start condition, the diffuser is open 
and is automatically closed immediately after the flow is established 
to provide maximum pressure recovery. 

The operation of the wind tunnel is automatic. All valves, gates 
and diffusers are hydraulic, and may be operated through a series of 
micro-switches and relay circuits by one button. By turning the 
selector switch to manual position, any of the operations can also be 
controlled individually, if desired. 
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Figure A-5 shows the details of the control console for the 
tunnel, along with an optional tape recorder. Notice the simplified 
controls and lighted control diagram mid -way of the center panel. In 
this visual-control schematic, each control element (valve, diffuser, 
etc.) is illuminated as it is activated either manually or electronically. 
The close proximity of the tunnel control console can be observed in 
Figure A- 1 . 

Additional information on the 7 - inch Bisonie Wind Tunnel can be 
obtained from reference 7. 


A- 2 Models 


This section describes the design characteristics of the models 
used in this test series. The model most frequently used in this series 
is a thin two-dimensional plate for generation of turbulent boundary 
layer and wake flows. The second model consists of a two-dimensional 
wedge for simulation of base recirculating flows. 

A-2.1 Thin Plate Model 


This model is essentially a two-dimensional turbulence- 
generating plate (Figure A- 6) constructed of aluminum. The leading 
edge has a 9°32 1 taper which extends back for 1.5 inches. The plate 
consists of a 1/4" thick flat surface for an additional 4.5 inches. 

The overall plate, 17 inches long, has a trailing edge of 11 inches 
with a very slow l°34 l taper. This model was mounted parallel to the 
floor on the horizontal centerline of the test section (Figure A-7) . 

It spans the full width of the tunnel for the first 5-1/4 inches, but 
has a 5/ 32- inch cutout or indentation on either side, avoiding full 
contact with the window of the tunnel (see Figure A-4) 4 The model show- 
ing surface finish, leading edge, and threaded holes (six 8-32) for 
mounting is photographed in Figure A-8. 


During these tests the models were operated with the Mach 
1.99 nozzle. The flow field resulting from this combination of nozzle 
and thin-plate model is depicted in Figure 2-3. 

The model was mounted as shown in Figures A-4 and A-5 . In 
this position, the sonic throat moves slightly downstream (positive x- 
direction) from the nozzle geometric throat, and dual sonic throats 
(one above and one below) are formed on the model about one to one and 
one-half inches downstream of the nozzle throat. It can be seen from 
the shadowgraph in Figure 2-3 that no shock rhombus is formed on by 
leading edge or reflection from the walls. 
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Assuming the throat location described above, one -dimensional 
Mach number, static pressure, density and temperature variations were 
calculated for the model-nozzle combination. This information is 
depicted in Figure A-9 which also shows a drawing of the geometric rela- 
tionship of the model and nozzle contour. The freestream flow speed in 
the test section is approximately 1660 fps. This information was used 
for comparison to convection speeds measured by the remote sensing 
system. 


A-2.2 Two-Dimensional Wedge Model 

The two-dimensional wedge model for simulation of recirculat- 
ing base flows is presented in Figure A-10. The model, about 7 inches 
wide, about 4 inches long, and about 0.75 inch thick, spans the full 
width of the tunnel. Constructed of aluminum, the model is relatively 
smooth and has a surface finish (about 125). The leading edge has a 
9° taper. The wedge is attached to the tunnel side walls by four 8-32 
screws with the axis mounted along the centerline, and the tapered 
leading edge facing upstream. The flow field generated by this model 
is shown in Figure 2-1. 


A- 3 Instrumentation 


This section describes the instrumentation used for sensing in this 
test series. 

A typical laser-schlieren wiring schematic for remote -sensing 
instrumentation is shown in Figure A-ll . Several variations of 
this arrangement are possible, depending on the type of flow problem 
and data desired (e. g. , amount of filtering, amplification, general flow 
conditions, etc.). The basic components as seen from the schematic are, 
in order, a source (or laser), photodetector, photodetector power 
supply, amplifiers, low and high cutoff filter (optional for some 
cases), oscilloscope for observation of raw data signal, analog corre- 
lator, and finally an oscilloscope and Polaroid camera for viewing and 
recording the final correlogram of the signals. 

Figure A-12 is a photograph showing part of the instrumentation 
installed at the 7 -inch Bisonic Tunnel. This photograph shows the 
general arrangement of the light source (placed for parallel beam 
measurements), the traversing and elevating stand, the analog correla- 
tor, and part of the conventional schlieren system from which flow studies 
can be made independent of the remote-sensing system. Figure A-13 shows 
more details of the electronic equipment used for processing the raw 
signals. This equipment is described in more depth in the following 
paragraphs . 
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A-3.1 Photodiodes 


The photodiodes used for signal detection were EG&G (Edgerton, 
Germmeshausen, & Grier, Inc.) Model S~D 100. This is a fast response 
light detector which operates in the visible and near-infrared spectrum. 
This silicon photodiode, operating over a wide spectral range, has fast 
response time combined with high sensitivity. Spectral response is from 
0.35 to 1.13 microns (near ultraviolet to near infrared). The sensi- 
tivity is 0.25 p.a/jLiw at 0.9p,. Typical values of rise and fall times 
are 4 x 10~ 9 and 15 x 10" 9 sec. Noise equivalent power (NEP) is in 
order of 1 x 10~ 12 watts/ \Tcps at 1000 cps. A dark current lower than 
0.2 x 10“ s amperes can be achieved with a bias of 10 volts. 

The temperature operating range varies from -65 to +100°C. 

The window is 0.12 M in diameter, constructed of Corning 7052 glass. 
Sensitive area is 0.11 in 2 (0.073 cm 2 ). Typical characteristics of 
this photodiode are summarized in the following table. 


TABLE A-l. PHOTODIODE CHARACTERISTICS 



Minimum 

Typical 

Maximum 

Units of Range 

Spectral Response 

0.35 


1.13 

Microns (10% points) 

Sensitivity (0.9p) 


0.25 


pa/ pw 

Rise Time 


4 x 10" 9 


Sec @ 90 

Fall Time 


15 x 10~ 9 


Sec @ 90 

Linearity 


R n 


n = 1.0 (within 5% 
over 7 decades) 

Bias 

1 


150 

Volts 

Capacitance 


8 


Picofarads @ 90 

Saturation Photo -Current 



0.12 

Peak amps from photo 
active area, load 
< 500 ohms at 90 bias 

Operating Temperature 

-65 

. i 


+ 100 

°c 
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During testing, these photodiodes are mounted in a convenient 
housing which has a BNC connector. The photodiode housing mounted on 
Uni-Slide elevating mechanisms are shown in Figure A- 14. The adhesive 
tape showing on the photodiode holds safety razor blades in position to 
split the beam, thus creating a quasi-schlieren effect (see section 
3.1.1). The photodetectors are mounted on an elevating and traversing 
stand similar to the laser source (Figure A-12) . Figure A- 15 shows the 
photodetector and mirror arrangement used for the cross beam measure- 
ment discussed in section 3.4.1. 


A-3.2 DC Power Supply 

The photodiode power supply consists of a 22-1/2 volt 
battery with an adjustable potentiometer to regulate the dc current. 
Normal operation is less than 1/2 unit output with a total power up 
to 1 milliwatt. Figure A-16 shows typical noise voltage as a function 
of frequency. These components are housed in a small metal box (2” x 
4” x 4") with appropriate BNC connectors attached for convenience of 
operation. 


Figure A-16 shows that, for a 22-1/2 V battery, this photo- 
diode has less than 5 percent noise from 500 Hz to large values of fre- 
quency. Figure A-13 shows the photodiode power supply n one line,” 
installed in its normal location. 

A-3.3 Amplifiers 


The amplifiers used for these tests were two (one per 
beam) Red cor Model 500, a two -stage type. The operating power of these 
amplifiers is 115V ± 10 percent, 60 cycle ac current at 10 watts . The 
gain for each stage ranges from a minimum of 10 to a maximum of 1000 in 
10 steps with accuracy of ±0.02 percent (dc) , linearity ±0.01 percent 
(dc) and stability ±0,01 percent (dc) . Bandwidth at 3 db full scale is 
1000 KC maximum (dc) and 10 cps minimum (dc) . Temperature operating 
range is 0 to 50°C. The following noise for the REDCOR is quoted with 
99.9 percent confidence. 


TABLE A-2 . TYPICAL AMPLIFIER NOISE 


Frequency 



Noise 


100 

KG 

10 

pv 

RMS 



1 

KC 

20 |li v peak 

to 

peak 

300 

cps 

14 

pv 

peak 

to 

peak 

30 

cps 

4 

pv 

peak 

to 

peak 
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The output is ±10 volts, ±10 milliamperes current for both channels, and 
the settling time for a full-scale step input is 20 microseconds to 0.01 
percent of final value at 100 KC bandwidth. Additional features of this 
amplifier is an adjustable bandwidth, high input impedance, and a band- 
width that is unaffected by gain change. It also features a solid state 
chopper . 


Figure A-13 shows a pair of Redcor amplifiers on line; the 
corresponding wiring diagram is shown in Figure A-ll. 

A- 3. 4 Fil ters 


An Allison variable filter (Model 2BR) was used on each 
leg of the remote sensing circuitry (see figures A-ll and A-13) . This 
filter features two separate networks: a high-pass double K section 
filter (0 to 20 KC) and a low-pass double K section filter (0 to 20 KC) . 
The low pass frequency was not used for this test; rather, the low pass 
was limited to 30 KHz by the Redcor low- pass frequency cutoff or to 
100 KHz when no low-pass filter was used. Each filter has 2 controls: 
an active band switch which changes the cutoff frequency in octave steps 
and a multiplier dial which tunes the cutoff frequency over one octave. 
The cutoff frequency is that which is attenuated approximately 3 db from 
the minimum insertion loss. The minimum band pass without additional 
insertion loss is 1/3 octave. The attenuation rate of these filters is 
about 30 db per octave. Each was pre- tuned to provide a minimum phase 
shift for signals over the frequency range anticipated in these tests. 
Distortion is about 1 percent for a 10-volt input, and at 1 volt, it is 
less than 0.10 percent. This instrument is recommended for measurements 
to 120 db below 1 volt. 

A- 3 . 5 Light Source 

The signal source was supplied by two types of lasers: a 

Spectra- Physics (S.P. ) Model 130 Gas laser and a Quantum Physics (Q.P.) 
Model LS32 laser. These lasers are illustrated in Figure A-17 in the 
position for parallel beam remote-sensing In the Bisonic Tunnel. Fig- 
ure A-18 shows the quantum physics elevated by a Lab-Jack in order to 
obtain cross-beam turbulence flow measurements in the wake of the thin 
plate. 


The operational and design characteristics of these lasers 
are examined in more detail in the following: 

° Spectra-Physics Laser . The Model 130 S.P. laser has a 
helium-neon gas filled plasma tube 27.5 cm long with 
an inside diameter of 2.5 mm. The ends of this tube 
are terminated with optical, s chi ieren- free , fused 
silica Brewster’s angle windows which result in a 
plane polarized output. Wavelength of this output 
is in the visible red circa 6328A. 
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The standard instrument has a hemispherical resonator 
with one spherical reflector (30 cm radius) and one 
planar reflector. The resonator output power is 0.75 
millivolts CW (minimum) from the spherical end. 

Resonator reflectors are also of optical quality, 
schlieren-free, fused silica. Reflectors are multi- 
layer dielectric coated for at least 99 percent 
reflectivity (at desired wave length) and are anti- 
reflection coated on back surfaces. 

Beam diameter is 1.4 mm at the exit aperture, and 
diverges less than 0.7 mil lir ad tans (145 seconds of 
arc). Power input requirements are provided by a 
self contained discharge exciter dc power supply. 

This requires a 60 cps 115 volt input of approximately 
90 va. The laser mode of operation normally used for 
these tests is TEM^. 

° Quantum Physics Laser . The Quantum-Physics gas -filled 
helium-neon laser (model LS32) used in these tests has 
a beam diameter of 1 mm, a power level output of 0.75 raw, 
and a wave length of 63281. In general this laser is 
similar to the Spectra-Phys ics Model 130 previously 
described. 


A-3.6 Miscellaneous Laboratory Equipment 

Other items of equipment consisted of two conventional 
Tecktronix Model 502A dual beam oscilloscopes, one for observing and 
measuring peak-to-peak voltage of the raw signal and one for viewing 
and recording the processed or correlated signal from the Princeton 
analog correlator. The latter was fitted with a swing-away Tektronix 
camera adapter. The signal was recorded on speed 3000 film by a 
laboratory Polaroid Land Camera. By using a manual shutter control, 
with an f setting of 16 and time exposure, the signal was traced on 
exposed film from a remote switch on the Princeton correlator where 
it had been temporarily stored in a multiple capacitor bank. A 
DIGITEX digital-type DC voltmeter was used to optimize power output 
by fine adjustment of the photodetector before each run. This equip- 
ment is shown in Figures A- 1 1 , A- 12, and A- 13. 


104 



A-4 Test Procedure 


Before each series of runs, the equipment was checked to insure 
that each instrument was adjusted properly. Also, the band pass was 
selected commensurate with the particular test objectives, and filters 
were phase-matched by use of a signal generator and oscilloscope. The 
next step in the procedure was to align the laser source in the desired 
geometrical position and focus the beam onto the photodiode until the 
photodetector power output was maximized. This power output was approxi- 
mately half of that available, since the remaining half of beam intensity 
was blocked by the knife-edge. 

Pre-run checks were made with and without flow to determine the 
extent of extraneous noise which was correlated in the signals. When 
ground loops and extraneous mechanical vibrations were minimized, the 
desired band pass was set, the gains of the amplifier were adjusted to 
avoid clipping the signal, and the correlator gains were adjusted to 
optimize the correctable signal to overcome correlator noise, yet at 
the same time avoid signal clipping. Blinking red lights warned when 
the correlator was overdriven. 

Direct current signals were recorded before each run, and the dc 
meter disconnected from the circuit to avoid ground loops. Atmospheric 
pressure, ambient temperature, and the inches of vacuum were also 
recorded before each run. The tunnel was activated by the automatic 
start button allowing about five seconds for the tunnel to settle to 
steady state operation. The visual display panel indicated the sequence 
of valves and outlet diffuser positioning. Visual observations of these 
were made during each run. Run time varied from about 40 to 60 seconds 
for each case. Data were correlated "on-line” and temporarily stored 
in the analog correlator. Upon completion of the run, the correlation 
function was displayed on the oscilloscope and recorded on a Polaroid 
photograph. 


A-5 Data Reduction Equipment 

Data taken during this test series were reduced on a Princeton 
Model 100 analog correlator. The primary function of this computer is 
to solve the following relationships: 


T 

lira tj; P fi (t) * f 2 (t - t) dt 

T^CO J 

where "C 11 is a cross-covariance of the fluctuating portion of random 
signals which contain flow information. Signals f x and f 2 , respectively. 
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represent information from two beams of electromagnetic radiation. If 
f 1 and f 2 are identical, then the result is an autocorrelation where the 
following result is obtained: 


C(t) = 


i • 1 

Inn - 

T-> 00 


f(t) • f(t - t) dt. 


This instrument operates as a hybrid computer to solve either of 
the two integrals for 100 n points of incremental time delay. At. 

The "nth" point is approximated by 

t-t 1 

RC 

C n^ ~ K / V A< t ’> ' V B (t ' - nt) &dt - 


where the following are true: 

° n RC n = time constant of the averaging circuit 
° "hAt" = time coordinates of the computed points 
0 t f represents past history (i.e . , t f > t). 

The three basic operations performed by the computer are (1) time 
shifting, (2) multiplication, and (3) integration. 

The computed values are stored in a 100-channel analog memory and 
may be recorded during and after computation on stripcharts, x-y 
recorder or an oscilloscope-camera combination. The latter method was 
used to record data for this report. Computed accuracy should not exceed 
1 percent deviation from the idealized function for any of the 100 points 
calculated . 

Other typical specifications of the Princeton Correlator are as 
follows : 

° Useful frequency range = Dc to about 250 KHz. 

° Averaging time constant (normally 20 sec) can be varied from 
0.1 to 400 seconds. 
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° Blinking red lights warn operator when the dynamic range of the 
input amplifiers is exceeded. 

° Computed values can be stored in the correlator with a decay 
that will not exceed 300 mv in 10 minutes at 25°C. 

° The compilation error is less than 1 percent of true value. 

° Calibration accuracy = ±2 percent. 

° Gain = 0.01 to 5.0 in 1 5 2, 5 sequence (with a variation of 
less than 0.1 percent /hour ) . 

° Zero drift = ±10 mv/hour. 

° Linearity I 1 percent. 


° Delay range = 100 microseconds to 10 seconds in 1, 2, 5 sequence. 

° Delay (time base t) accuracy = ±1 percent at mid-range and ±2 
percent at five fastest and three slowest ranges. 
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FIG. A -6. DRAWING OF THE THIN PLATE MODEL 
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APPENDIX B 


COMMENTS ON THE EXPERIMENTAL NOISE STUDY IN 
MSFC'S 14-INCH TRISONIC WIND TUNNEL (TEST # TWT 395) 


The experimental investigation presented in this report originated 
with the objective of conducting remote-sensing cross-beam experiments 
in MSFC’s 14 x 14-Inch Trisonic Wind Tunnel to demonstrate the capabil- 
ity of measuring supersonic turbulence parameters using a laser source. 
The Trisonic Wind Tunnel cross-beam tests were conducted in the 14-inch 
tunnel special test section which is basically an annular nozzle with 
provisions for attaching an axisymmetric blunt-base body or a nozzle 
and vertical and horizontal windows for cross-beam access. Subsequent 
analysis of these cross-beam data showed essentially negative results 
because of an unsteady flow field in the Trisonic facility. 

The major problem of the Trisonic facility was the almost complete 
domination of the data by facility- induced noise. This appendix was 
written to support the contention that excessive facility-induced mech- 
anical and acoustic noise dominated the data and thus precluded fruit- 
ful cross-beam tests in the facility (at least until the remote-sensing 
system has been further developed). The reasons why the 14- inch tunnel 
was not satisfactory for remotely sensed measurements are discussed in 
the following paragraphs. 

Figure B-l shows fifty-seven qualitative correlations of various 
beams and accelerometers which were located both inside and outside the 
turbulent flow field of interest. The first three rows across the top 
depict the autocorrelation of the individual horizontal and vertical 
beams and horizontally and vertically oriented accelerometers which 
were mounted on the outer wall of the STS. Sixty percent of these show 
definite periodic trends which are detrimental to the extraction of use- 
ful turbulent flow information. 

The four remaining rows across the bottom depict cross-correlation 
of the horizontal and vertical beams and cross-correlation of the single 
beam signals with the signal from the individual accelerometers. This 
latter device gives a qualitative indication of the relative magnitude 
and frequency of noise common between the flow field and structure. As 
for the autocorrelations, these cross-correlations are dominated by an 
extraneous periodic noise. Attempts to filter these data were not 
successful because of an adverse signal- to-noise ratio. 
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Figure B-l, Column (1), shows an attempt to obtain an intensity 
traverse. This should have resulted in a correlation function with a 
minor peak at zero time l^g since the beams are crossed in a region of 
low turbulence. However, the peak was masded by a periodic signal of 
1380 cps. A great deal of this noise can be attributed to structural 
vibration. This is indicated by the cross -correlation of the accelero- 
meter with each individual beam, since definite periodic correlations of 
significant magnitudes apparently exist between them. Also it may be 
noticed that these beams were going through a plastic sleeve extension 
of the nozzle outer wall which probably contributed to the correlation. 
Moving the beams downstream two inches (out of the plastic sleeve) 
apparently did not improve the situation significantly (see Column 2). 
The flow conditions and coordinates of column 3 are the same as (1) and 
(2) except a dominant peak should have been obtained at zero time lag, 
since the correlation volume is located in a dense shear layer. In 
Column (4) the correlation showed a definite mechanical or acoustical 
vibration between beams because one beam was completely removed from 
the flow and directed along side of the tunnel wall. Furthermore, 
Columns 4, 6, 7 and 8 depict various attempts to soft-mount the lasers 
and stands, but this apparently did not eliminate the vibrations. 

Upon analysis of these data, the following conclusions and recom- 
mendations were made: 

° The flow field, structure and surrounding environment in the 
14-inch TOT (STS) facility were dominated by an extraneous 
periodic noise which precluded ascertaining useful data at 
the time of the test. 

° Attempts at filtering these raw data signals and various 
methods of soft-mounting the lasers and detectors did not 
alleviate the situation so that flow related correlation 
could be obtained. 

0 Remote sensing measurments in the 14 x 14- inch wind tunnel 
might be possible after (1) more analysis of the 14-inch 
TOT noise problem, and (2) development of techniques to 
improve signal- to-noise ratio of the remote sensing system. 

° A different facility was needed with a lower facility- induced 
noise level. 

These results essentially led to the selection of the 7 x 7-inch 
BWT facility for the immediate future development of the remote-sensing 
tool. This approach appeared to offer the most economical method of 
expediting development of the remote-sensing system for measuring super- 
sonic turbulence. 
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APPENDIX G 


UNKNOWN ANOMALIES IN DATA ACQUISITION SYSTEM 


Unknown anomalies should be expected in any research work; the 
research involves identifying and resolving them. The anomalies in 
the system were (1) low-pass filters set too low, thus blocking signi- 
ficant amounts of data; (2) use of attenuation after amplification, 
which, combined with low maximum amplifier input voltage, probably 
resulted in clipping of signals; (3) high-pass filter set too low, 
thus allowing undesirable flow noise and facility- induced noise in 
the signals; (4) filtering low frequencies (e.g. , 60 eye/ sec) with 
large amplitudes before the run, probably causing clipping of signals 
during the run, which was not detected after the clipped signal was 
filtered; and (5) photodetectors and other possible points of contact 
between grounds were not isolated (not always a source of trouble). 

As can be seen, the problem areas were associated with the ampli- 
fiers (i.e. , clipping) and filtering methods . 
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